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1.0 SUMMARY 


A ROTO architecture, braking and steering control law and display designs for a research high 
speed Rollout and Turnoff (ROTO) system applicable to transport class aircraft are 
described herein. Minimum surface friction and FMS database requirements are also 
documented. The control law designs were developed with the aid of a non-real time 
simulation program incorporating airframe and gear dynamics as well as steering and braking 
guidance algorithms. 

An attainable objective of this ROTO system, as seen from the results of this study, is to 
assure that the studied aircraft can land with runway occupancy times less than 53 seconds. 
Runway occupancy time is measured from the time the aircraft crosses the runway threshold 
until its wing tip clears the near side of the runway. Turnoff ground speeds of 70 knots onto 
30 degree exits are allowed with dry and wet surface conditions. 

Simulation time history and statistical data are documented herein. Parameters which were 
treated as variables in the simulation study include aircraft touchdown weight/speed/location, 
aircraft CG, runway friction, sensor noise and winds. After further design and development 
of the ROTO control system beyond the system developed in reference 1, aft CG MD-1 1 
aircraft no longer require auto-asymmetric braking (steering) and fly-by-wire nose gear 
steering. However, the auto ROTO nose gear hysteresis must be less than 2 degrees. 

The 2 sigma dispersion certified for MD-1 1 CATIIIB is acceptable. Using this longitudinal 
dispersion, three ROTO exits are recommended at 3300, 4950 and 6750 feet past the runway 
threshold. The 3300 foot exit is required for MD-81 class aircraft. 

Designs documented in this report are valid for the assumptions/models used in this 
simulation. It is believed that the results will apply to the general class of transport aircraft; 
however further effort is required to validate this assumption for the general case. 
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2.0 INTRODUCTION 


OBJECTIVES 


The objective of this study was to design a research ROTO guidance and control system 
which will be used to develop operational ROTO requirements for both automatic and 
manual piloted operation under normal and reduced visibility conditions. This ROTO 
system will utilize the satellite navigation system called Global Position System (GPS) 
operated in a differential mode. 


The study addresses the design of : 

a variable auto-braking control law, 
a steering control law, 

ROTO exit prediction logic, 

ROTO arming procedures, and 
ROTO displays. 


Another objective is to develop the runway/exit friction and navigational database 
requirements. 

Manual ROTO steering and deceleration guidance control laws were not designed as part of 
this study. 

REQUIREMENTS 


The requirements used in this study are outlined in reference 1 . The requirements were: 

Nose Wheel/Rudder Steering 
Braking/Reverse Thrust 
Touchdown Dispersion 
Navigation Noise and Update Rates 
Longitudinal/Lateral Acceleration and Jerk 
Control Law Update Rates 
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3.0 ROTO CATniB ARCHITECTURE 


A preliminary assessment of the impact of CATIIIB requirements on the High Speed ROTO 
system requirements has been made and is herein described. The first step was to create a 
preliminary Functional Hazard Analysis to assign levels of criticality to the various functions 
and components of the ROTO system. A top level system diagram was prepared to 
facilitate the analysis and identify the main hardware or software elements of the system. 

The results of the FHA were used to help determine the level of redundancy of various 
system elements and identify special monitoring or procedural requirements. The results of 
the following Functional Hazard Analysis suggest that most of the ROTO system 
components supply a critical function and erroneous outputs or loss of function at a critical 
time can produce hazardous or even catastrophic results. In order to satisfy FAA and other 
safety requirements it is therefore necessary to show that no single failure or combination of 
failures not shown to be extremely improbable will produce these results. In some cases this 
is accomplished with simple redundancy and monitoring. In other cases where adequate 
redundancy alone is impractical, it will be necessary to rely on showing that a more modest 
level of fault tolerance combined with a short exposure time to the critical case is sufficient. 

In fact for much of the system components and interfaces, the redundancy and monitoring 
required to satisfy the concerns elicited by the attached FHA are already provided by the 
basic CATIIIB autoland system. Only the new components/tasks imposed by auto and 
manual High Speed ROTO function are addressed here. The shaded areas in figures 3.1 and 
3.2 identify these items which are discussed below. 

It is assumed that any software/hardware failure on the runway would revert to non-ROTO 
and a failure on the exit would remain in its current mode; auto or manual. The pilot would be 
free to take over from auto ROTO if he chooses. 


FMS Data Base 


The accuracy of the airport runway and taxiway geometric information provided by the FMS 
is critical to the ROTO guidance commands. In addition to assuring the data accuracy when 
loading the data into the database, detection of any failure affecting the data during use will be 
necessary. It would be unrealistic to assume more than two FMS functions can be provided. 
If there are only two FMS data bases, it is not clear how this can be construed fail operational 
unless all necessary data is transmitted to the ROTO computers prior to the start of the 
autoland after being validated by cross checking. FMS data may also contribute to estimating 
aircraft weight, drag and thrust characteristics; which are inputs to the ROTO software exit 
prediction logic. 

[A dual FMS and data base is assumed.] 
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GPS 


Positioning data provided by the DGPS receiver and associated local area augmentation signal 
source must be fail operational and provide the necessary reliability level to meet the 
"extremely improbable" criteria for loss of function. It is assumed this would also be the 
source of position data for the autoland guidance. No additional integrity requirements would 
be imposed as a result of the ROTO function. 

[The DGPS airborne elements are assumed to be triplex with outputs voted within 
each using computer (ROTO and/or FCC).] 


VHF Communications Link with Ground Station 

Transmittal of estimated runway friction coefficient to the approaching aircraft is not a 
critical element since its absence would not affect the safety of the rollout. However, its 
accuracy and integrity must satisfy some minimal reliability requirement in order to take 
credit for advance prediction of the turnoff to be used. 

The integrity of the voice communications with the control tower may be more critical since 
it is conceivable that a last minute abort may be required as a result of a taxiway 
encroachment incident. The tower would be expected to detect this situation and alert the 
approaching aircraft. The combination of the probability of this event and the probability of 
loss of VHF communications should be sufficiently remote that no new requirements are 
imposed on the communications link. 

[VHF communications are assumed to be dual.] 


ROTO Computer 


All steering and deceleration commands originate in this unit both for automatic and manual 
mode of operation. The accuracy of commands and continuity of function must be assured 
throughout the rollout. This function would likely be incorporated within the same 
computers as the auto land functions (viz.; flight control computer), although it has been 
shown as a separate element on the reference diagram for clarity of function. 

[A dual-dual or triplex implementation is assumed.] 

Servo Backdrive for Rudder Pedals and Tiller 


For pilot situational awareness, backdrive of all control functions except braking has been 
deemed to be necessary during the automatic ROTO mode. For the rudder pedals this should 
impose no new requirements since it would also be required for normal CATIIIB autoland. 

In addition to backdriving the tiller for situational awareness, auto ROTO now requires the 
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ROTO computer to command the tiller via a dual servo which is fail passive (similar to MD- 
80 autopilot servo). The tiller signal compensates for the hysteresis in the rudder to nose 
gear mechanical control system. A dual fail passive servo is deemed adequate because a servo 
failure on the runway would revert to non-ROTO. Unlike other failures, this particular tiller 
servo failure on the exit would revert to manual ROTO, due to the undesirable MD-1 1 auto 
ROTO nose gear hysteresis greater than 2 degrees without tiller compensation. 

[Dual backdrive for rudder pedals and simplex for tiller backdrive is assumed.] 


Servo Drive for Reverse Thrust Control 


If reverse thrust is judged to be necessary for adequate deceleration performance in automatic 
mode, then the position of the throttle levers will have to be moved to reflect the engine 
commands even with reversers deployed. This could be a significant complication unless a 
clever mechanical design of the throttle quadrant is contrived that allows the autothrottle 
servo[s] to move the levers while in reverse thrust command position as well as normal 
position. 

[A dual throttle servo is assumed.] 


Braking System 


Braking action in response to variable automatic or manual commands must be dependable 
and consistent in order to achieve the deceleration profiles required to make the designated 
high speed turnoffs. Some current aircraft have dual braking systems. This should be 
sufficient to satisfy reliability and safety requirements if it can be shown that there is no loss 
of control for a single system failure and that the system adjusts adequately to the reduced 
braking pressure available. The braking computer is also required to be fail operational. 

[A dual braking system is assumed.] 


Nose Gear Steering Actuation 


This (along with baseline rudder control) is the most critical function of the high speed 
ROTO system. Failure of steering function at a critical moment would probably be 
catastrophic. All portions of the control path from and including the ROTO computer/FCC 
to the actuators and connecting linkages must be shown to be fail operational. An anti- 
hysteresis software/tiller cable grabber system will reduce MD-1 1 auto-ROTO nose gear 
hysteresis down to an acceptable 2 degrees. 

[A dual installation is assumed.] 
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ROTO Display(s) 


If the ROTO maneuver is required to be initiated manually, then the every element required 
for satisfactory control must be fail operational. This includes the mechanism for display of 
runway map and/or command guidance. Loss of displayed information at a critical time is 
unacceptable. Dual display mechanisms are certainly a requirement. Furthermore, it must be 
shown either that the loss of the pilot's display at a critical moment is extremely improbable, 
or that it is possible to transfer control between pilots without significant disturbance even 
during the critical time. 

[A dual display arrangement with fail obvious characteristics is assumed]. 


A preliminary reliability analysis has been performed to show that the redundancy levels 
selected for critical system components is adequate. This analysis is included following the 
FHA. In addition to the simple component redundancy and theoretical reliability analysis, it 
will also be necessary to show that no single point failures or common failure modes exist 
that could invalidate the redundancy. For example, possible software errors must be 
eliminated either by stringent design and test measures or by application of dissimilar design 
techniques or a combination of both. 
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PRELIMINARY FUNCTIONAL HAZARD ANALYSIS 


INTRODUCTION 


This preliminary functional hazard analysis (FHA) was prepared in accordance with DAC 
Design Safety Manual PD-503 and FAA Advisory Circular 25. 1309-1 A. It is an initial risk 
assessment of the functions and evaluates and identifies critical safety areas for potential 
hazards. This information is used in developing system safety requirements and establishing 
the framework for other analyses. While the FHA is not a compliance document, it is the 
basis for certain design decisions and does provide guidance for decisions on other analyses 
for certification. 

The purpose of the FHA is to develop safety design requirements for the system(s) which 
will perform the function being considered and establish the framework for the Certification 
Plan and subsequent assessments. It provides information about potential functional failure 
conditions which should be used for establishing the required system architecture, software 
integrity level requirement, system separation and isolation requirement, and minimum 
equipment list (MEL) requirements. The subsequent design and development process must 
result in a system which, as installed, can be shown to meet applicable FAR Part 25 safety 
standards. 

SYSTEM DESCRIPTION 


The baseline high-speed rollout and turnoff (ROTO) guidance and control system architecture 
approach is shown in figure 3. 1 and 3.2. The ROTO system architecture provides for an 
automatic ROTO mode as well as a manual ROTO mode. Both modes require the following 
equipment/functions : 

1 . A ROTO navigation input which is provided primarily by a global positioning system 
(GPS) operated in a differential mode and supplemented by an inertial reference 
system. 

2. Ground operational data such as available exits, exit speeds, winds, temperature, 
surface conditions, other aircraft, runway friction, aircraft weight, etc. 

3. ROTO exit locations/geometry which is permanently stored in the flight management 
system (FMS). 

4. Control panel for the selection or de-selection of ROTO automatic or manual mode as 
well as ROTO abort. 

5. ROTO displays that provides the ROTO guidance information on the head-up- 
display (HUD) and/or primary flight display (PFD). 
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6. Flight control computer to provide nose gear steering control, rudder actuation, and 
feedback signals to backdrive the tiller, rudder pedals, and throttle levers. 

7. Brake computer to provide the braking function 

8. Reverse thrust function. 

9. ROTO computer 

In the automatic ROTO mode the following operational sequence of events occurs: 

1 Air traffic controller (ATC) offers all available ROTO exits to the ROTO software. 

2. ROTO software offers the pilot one ROTO exit. 

3. Pilot arms (accepts) ROTO exit. 

4. ROTO deceleration provided by the auto braking and auto reverse thrust 

5. Auto ROTO steering using the rudder, and nosewheel steering. 

6. Pilot monitors the ROTO automatic mode by the HUD and/or PFD. 

In the manual mode, the following operational sequence of events occurs: 

1 . ATC offers all available ROTO exits to the ROTO software. 

2. ROTO software offers the pilot one ROTO exit. 

3. Pilot arms (accepts) ROTO exit. 

4. Pilot brakes, commands reverse thrust, and steers onto the ROTO exit by following 
the ROTO guidance information provided by the ROTO displays either on the HUD 
and/or PFD. If the pilot preselects auto brake/auto reverse thrust manual ROTO 
mode, the sequence 4 event is changed and an event 5 is added as follows: 

4. ROTO deceleration provided by auto braking and auto reverse thrust 

5 . Pilot steers onto ROTO exit by following the ROTO guidance information provided 
by the ROTO display either on the HUD and/or PFD 


Refer to figure 7.2 for a pictorial representation of the automatic and manual ROTO modes 
operational concept. 

Preliminary ROTO guidance and control system architecture definition has identified the 
ROTO computer, the backdrive servos for the tiller, throttle levers in the reverse thrust 
mode, and variable auto-braking as the major unique ROTO system hardware additions to the 
baseline large and/or heavy transport aircraft. Modifications to the baseline aircraft system 
equipment are anticipated to be needed also to meet reliability and safety requirements, but 
the specific modifications have not been determined at this time. 
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SYSTEMS FUNCTIONS 


The initial step in the development of the preliminary ROTO FHA consists of the 
identification of the ROTO system functions. The following functions list has been 
established for the preliminary FHA: 

1. ROTO navigation 

2. Runway acquisition clearance 

3. ROTO steering guidance 

4. ROTO deceleration 

5. ROTO displays 

HAZARD CLASSIFICATION 


The hazard classifications used in this FHA are derived from requirements set forth in the 
Code Of Federal Regulations (CFR) and are defined as follows: 


Class IV (Minor): Failure conditions which would not significantly reduce airplane safety, 
and which involve crew actions that are well within their capabilities. Minor failure 
conditions may include, for example, a slight reduction in safety margins or functional 
capabilities, a slight increase in crew workload, such as routine flight plan changes, or 
inconvenience to occupants. 

Class III (Major): Failure conditions that would reduce the capability of the airplane or 
the ability of the crew to cope with adverse operating conditions to the extent that there 
would be, for example, a significant reduction is safety margins or functional capabilities, 
a significant increase in crew workload or in conditions impairing crew efficiency, or some 
discomfort to occupants. 

Class II (Hazardous): Failure condition which would reduce the capability of the airplane 
or the ability of the crew to cope with adverse operating conditions to the extent that 
there would be, for example, a large reduction in safety margins or functional capabilities, 
higher workload or physical distress such that the crew could not be relied upon to 
perform its tasks accurately or completely, or adverse effects on occupants. 


(Note: AC 25. 1309-1 A refers to this condition as "severe major.") 

Class I (Catastrophic): Failure conditions that would prevent continued safe flight and 
landing. 
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AC 25.1309-1 A defines continued safe flight and landing as: 


"The capability for continued controlled flight and landing at a suitable airport, possibly 
using emergency procedures, but not requiring exceptional pilot skill or strength. Some 
airplane damage may be associated with a failure condition during flight or landing" 


SUMMARY OF RESULTS/RECOMMENDATIONS 


The following summary of results/recommendations are based upon the functional hazard 
analysis worksheets contained in Appendix A beginning on page A- 133. Tables 3.1 and 3.2 
present the summary of Class I and Class II hazards. The FHA worksheets have been 
prepared based upon the high level ROTO guidance and control system functional and 
operational description presented in the SYSTEM DESCRIPTION section above. The 
objective is to identify top level critical system hazards and provide a high level system 
design approach to address the identified hazards in terms of system architecture 
recommendations/ considerations, such as selective application of redundancy, increased 
equipment reliability, failure monitoring concepts/schemes, etc.. Based upon tables 3.1 and 
3.2 summary of Class I and Class II hazards, the following system architecture 
recommendation/.considerations have been developed and are provided in tables 3.3 and 3.4. 
Tables 3.3 and 3.4 include the system function, hazard description, baseline aircraft or unique 
ROTO system critical equipment that contributes to the Class I or Class II hazards, and the 
specific system architecture recommendations/ considerations for each identified critical 
system. 


CONTINUING STUDIES/ANALYSES 


One of the primary tasks of the reliability and safety effort for the design development 
process of the ROTO guidance and control system architecture will be the continuing 
development of the quantitative and qualitative tradeoff parameters for the selection process 
of fault-tolerant architecture. This initial preliminary FHA and the summary of system 
architecture recommendation/considerations for Class I and Class II hazards presented in 
tables 3.3 and 3.4, form the basis for these continuing studies/analyses. As mentioned in the 
preceding System Description section, modification to the baseline aircraft system equipment 
have not been determined at this time, but the development of theses unique modifications 
will be part of the continuing studies/analyses effort. In addition, the continuing reliability 
and safety studies/analyses will provide the reliability and safety design criteria for 
establishing the fault tolerant system architecture, software integrity level requirement, 
system separation and isolation requirement, and minimum equipment list (MEL) 
requirements that are part of the process in meeting the system safety certification 
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requirements contained in the Federal Aviation Administration (FAA) Code of Federal 
Regulation (CFR) Chapter I, Federal Aviation Regulation (FAR) Part 25.1309. 

Examples of some of the ROTO guidance and control system architecture and critical 
equipment issues that may be addressed in future studies/analyses include: 

1 . Selective application of reliability and safety enhancements such as redundancy, 
increased equipment reliability, etc. for the following: 

a. Global positioning system 

b. ROTO computer 

c. VHF data and communication channels 

d. Flight control computer 

e. Brake computer 

f. ROTO display electronics that drives the HUD and/or PFD 

2. Failure monitoring concepts/schemes that addresses: 

a. Erroneous navigation data 

b. Erroneous deceleration commands 

c. Erroneous reverse thrust 

d. Erroneous steering guidance 

e. Erroneous rudder 

f. Erroneous nosewheel steering 

g. Erroneous aircraft weight information 
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ROTO SYSTEM RELIABILITY ANALYSIS 


A preliminary reliability analysis was prepared based upon the current high speed rollout and 
turnoff (ROTO) guidance and control system description. The preliminary ROTO system 
functional hazard analysis established that the ROTO system quantitative safety certification 
requirement shall meet the extremely improbable case as defined in the Federal Aviation 
Regulation (FAR) Part 25.1309 and the Advisory Circular 25.1309.1A as follows: 


The occurrence of any failure condition which would prevent the continued safe flight and 
landing of the airplane is extremely improbable. Extremely improbable failure conditions 
are those having a probability on the order of 1 x 10E-09 or less. 


This preliminary reliability analysis is an integral part of the overall ROTO system design 
development process by establishing the quantitative reliability probability estimate to 
validate that the proposed ROTO system architecture has the potential to meet the 
extremely improbable quantitative requirement. 

The results of the ROTO system reliability analysis indicates that the estimated probability 
of failure is 2.06 x 10E-10 as shown in table 3.5 ROTO SYSTEM PROBABILITY OF 
FAILURE. 

The reliability analysis was developed based upon the mean-time-between-failure values for 
the applicable equipment as shown in the table 3.5, the ROTO system reliability block 
diagram shown in figure 3.3, and an operational time of five minutes. 
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4.0 MODELING 


The model used in this study was outlined in reference 1 . The model was implemented in 
both FORTRAN 77 and MATLAB SIMULINK diagrams, which have been delivered to 
NASA Langley. The figure found in the appendix on page A-63 illustrates the MATLAB 
SIMULINK ROTO top level diagram. The aircraft simulation is a 3 degree of freedom (yaw, 
forward, lateral). It calculates aerodynamic, thrust and tire forces on the airplane and solves 
the resulting equations of motion to determine aircraft accelerations, velocities and positions 
during a simulated turnoff. The simulation also includes hydraulic models of the nosewheel 
steering, rudder and autobrakes. The simulation begins at touchdown. The model includes 
the following items: 


ROTO Exit Geometry 

Nosewheel, Rudder and Autobrake Actuation 
Tire-runway Coefficient of Friction 

Forces - Aerodynamic, Thrust, Braking Drag, Main & Nose Gear (Vertical & Side) 

Aircraft Equations of Motion - Acceleration, Velocity, Position 

Navigation 

Winds 

ROTO Control Laws 
Exit Prediction Logic 


An addition was made to the model to include rudder cable hysteresis (automatic ROTO) and 
nose gear cable hysteresis. The figure in Appendix A on page A-l 12 illustrates the cable 
hysteresis. The anti-hysteresis block is a software/hardware feature to minimize the auto 
ROTO nose gear hysteresis to 2 degrees. This feature is described in the Steering Control 
Law ROTO Design section 7. 
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A variety of aircraft types may be simulated by providing the simulation with unique aircraft 
characteristics (i.e. dimensions, actuator, gear, aero and thrust characteristics). Some of these 
characteristics are described below for an MD- 1 1 and MD-8 1 : 


Unique Aircraft Characteristics 

MD-11 

MD-81 

Landing Weight (1000 lb.) 

340 to 480 

82 to 128 

CG 

12% to 34% 

-0.8% to 33.4% 

Landing Air Speed (Keas) 

130 to 166 

110 to 143 

Wing Span (ft) 

165 

108 

Aircraft Yaw Inertia (million slug-ft A 2) 

26 

4.1 

Wing Gear: 

Number 

2 

2 

Wheels/Gear 

4 

2 

Center Gear: 

Number 

1 

0 

Wheels/Gear 

2 

0 

Nose Gear: 

Number 

1 

1 

Wheels/Gear 

2 

2 

Distance between wing and nose gear (ft) 

81 

72 

Distance between wing gear (ft) 

35 

17 
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5.0 DATABASE REQUIREMENTS 


It is assumed that long term permanent data (ROTO exit geometry data) would be stored in the 
Flight Management System (FMS) and temporary uplinked data (used in approach and roll-out) 
would be stored in the Flight Control Computer (FCC). 


FMS ROTO EXIT GEOMETRY REQUIREMENTS 


The MD-1 1 FMS database currently stores its permanent navigation data in 1 MB. Storage capacity 
is therefor at a premium. Each 3 ROTO exit airport would require approx. 400 bytes of FMS 
permanent storage. ROTO exit prediction logic requires approximately 20 aircraft parameters (4 
bytes each) which could be permanently stored in the FMS. 

1. The FMS database currently stores the runway threshold and end points in latitude/longitude 
coordinates. ROTO requires that the resolution of these coordinates be an inch. This would 
require the use of lOOOths of a second in polar notation or 7-8 decimal places using a 4 byte real 
number notation. 

2. ROTO requires that aircraft DGPS position accuracy be at least +/- 1 feet. The resolution of the 
aircraft location should also be an inch. If the resolution is as much as 1 foot, the accuracy should 
tighten to not be greater than +/- 1 foot. 

3. ROTO, developed in this study, uses an X,Y coordinate system to represent the runway and its 
exits. The runway is represented by the X axis with the threshold at the origin and the runway 
extending along the positive X axis. The right side of the runway, as it appears to a landing 
aircraft, has positive Y values and, therefore, a right hand exit would have positive Y values. The 
runway centerline has zero Y values. The units (resolution) of X and Y should be 1/2 foot. The 
accuracy of the X,Y coordinates should be an inch. The endpoint of a 12,000 foot runway would 
have coordinates of X = 24000 and Y = 0. The resolution of 1/2 foot allows X and Y to be 
represented by 2 byte integers, thereby minimizing X,Y pair storage requirements. 

4. An array of 2 byte integer X,Y coordinate pairs in units of 1/2 foot represents an exit as a series 
of line segments. It is assumed that the exits cannot necessarily be represented by exact formulae 
such as an arc. The length of any particular line segment must not cause any point of the line 
segment to have a perpendicular error from the actual exit path greater than 1/2 foot. Using this 
criteria requires about 22 X,Y pairs to represent a 2200 foot long ROTO spiral (reference 2) 30 
degree exit, reference 2. The exit centerline radius (2 byte integer, required by steering control 
law) is stored with each X,Y pair. Each 3 ROTO exit airport would require about 400 bytes of 
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FMS storage (3 * 2 bytes (X,Y,R) * 22 (2200 feet/- 100 foot spacing) * 3 (3 exits/airport)). The 
exit centerline radius values may be stored in units of feet. 


5. The aircraft s latitude/longitude position must be projected onto the runway's X,Y coordinate 
system using the following steps: (all trig functions must be accurate to 8 decimal places) 


a) Convert latitude/longitude positions from the FMS data base (polar coordinates) to 
floating point numbers. 



THRESHOLD 

LATt/LONGt 


DEG2FT = - * ^Omin 

min deg 

b) Calculate the angles e and a. 


Angle e only needs to be calculated once as follows: 

HE = abs (LONGt - LONGe ) *DEG2FT* cos|^— ATt * LATe ) j ( feet ) 


VE = abs (LATt - LATe ) *DEG2FT(feet) 


~\( HE 
e = tan 1 — 
VE 


Angle a would be calculated as: 


HA = abs (LONGt - LONGa )*DEG2FT* cosf ~ T[ + LATa ) V feet ) 
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VA = abs (LATt - LATa ) *DEG2FT(feet) 


HA 
VA 

c) Account for the affect of runway direction in defining left and right side of the runway. 

if (LONGe is to the east of LONGt) then 
sign = 1 

else 

sign = - 1 

endif 



d) Calculate the distance A from the runway threshold to the aircraft location. 
A = -VVA * *2 + HA * *2(feet) 


e) Calculate the aircraft's X and Y position (4 byte real number) in units of 1/2 feet. 


Aircraft's X position = A * cos(a - e) * 2 (1/2 feet) 


Aircraft's Y position = sign* A * sin(a - e) * 2 (1/2 feet) 


f) After calculating the aircraft’s lateral deviation/deviation rate from the runway/exit 
centerline, the lateral deviation/deviation rate should be divided by 2 to convert them to 
units of feet for use by ROTO control law algorithms. 
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FCC ROTO UPLINKED DATA REQUIREMENTS 


The uplinked data required by ROTO are DGPS ground-computed corrections for the continuous 
onboard update of the aircraft position and ground velocity and other data specific to the current 
approach (runway friction measurements, winds, available exits and exit speeds) for use by the 
ROTO exit prediction logic. If exit prediction were performed by ATC on the ground, the aircraft 
would only need to know its designated ROTO exit and exit speed. However, the aircraft weight 
may then have to be sent to the ground. All other data would still be required by the ATC. Memory 
space for approximately 20 aircraft parameters (4 bytes each) need to be reserved in the FCC 
(calculated) or FMS (permanent). A block diagram of the exit prediction algorithm is shown in figure 


Uplinked 

1 . All currently available ROTO exit locations and their accompanying exit entrance speeds. Exit 
geometries are permanently stored in the FMS navigation database. 

2. Friction. An array of speed dependent surface conditions at locations along the runway (approx. 
250 feet spacing). The surface conditions should be an average or represent the entire 250 foot 
length, otherwise exit prediction may result in false conclusions. The measurement spacing may 
not be as important as its accurate representation of the distance it covers. 

3. Uplinked estimated winds and onboard DGPS position and velocity data would be used to 
predict the aircraft touchdown position and ground speed. Estimated flare distance may depend 
on approach mode (autoland, ILS manual, non-ILS manual). The touchdown location of a non- 
ILS manual approach may be difficult to predict. Use of ROTO exit prediction logic in ROTO 
operation may require the engagement of ILS manual or autoland before ROTO can be engaged. 
Estimated crosswind is also used within the exit prediction logic. 


Onboard 


1. Known (modeled) Aircraft Parameters: 

Thrust: Arrays of the following thrust profiles: 
Approach Idle versus Time 
Approach Idle versus Airspeed 
Ground Idle versus Airspeed 
Reverse Idle versus Airspeed 
Maximum Reverse versus Airspeed 
Aircraft Wing Area 
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Percent MAC 

Distance from nose gear to CG 
Distance from main center gear to CG 
Distance from main wing gear to CG 
Aircraft CG Height 
Forward most CG 
Aft most CG 

Drag Coefficient at forward CG (based on aircraft configuration, surfaces deployed) 
Drag Coefficient at aft CG (based on aircraft configuration, surfaces deployed) 
Average Anti-skid Efficiency 

Expected Nose Gear T.D. Time (time starts after main gear touchdown) 

Reverse Thrust Available Flag 
Automatic Reverse Thrust Available Flag 
Maximum allowed deceleration 
2. Estimated Aircraft Parameters: 

Aircraft Weight 
Aircraft CG 
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6.0 FRICTION REQUIREMENTS 


This section will describe friction requirements which were found through simulation and a 
literature search. 


SIMULATION RESULTS 


The friction requirements from simulation results were developed with non-grooved runway 
and exit models. Grooved runways and exits have a higher friction coefficient than non- 
grooved ones and should be considered for the ROTO operation. 


Minimum Allowed Runway Friction 


The runway occupancy times documented in this report were determined with exit locations 
at 3300, 4950, 6750 and 8000 feet past the runway threshold, using the non-grooved dry and 
wet friction curves found in figure 6.4. Different exit locations may be required for lower 
surface friction conditions (e.g. flooding, and ice). The exit locations used in this study will 
accommodate some amount of flooding and ice less than 10% and 5%, respectively, of the 
wet friction condition. 

Some low-probability-of-occurrence, heavy/late MD-1 1 landings on a wet surface required a 
ROTO exit past 8000 feet and, thus, were considered for this study to be non-ROTO 
landings. 


70 Knot ROTO Exit 


For acceptable deceleration and steering on the 70 knot exit, the exit must have friction at or 
better than the following runway condition combinations: (see Maximum Ground Coefficient 
Friction for MD-1 1 Diagram, figure 6.4) 


90% Textured Concrete Wet & 10% Textured Concrete Flooded 

OR 

95% Textured Concrete Wet & 5% Ice (0.02*Dry) 

OR 

At or above an average constant friction coefficient of 0.2 (i.e. alternating 0.35 and 
0.05 patches smaller than 15 feet in length). The Minimum Allowed Runway Friction 
section above allows a friction coefficient less than 0.2 on the runway. The Textured 
Concrete Wet condition has a lower friction coefficient than 0.2 at runway speeds. 


23 



If we assume that the entire exit has the Textured Concrete Wet surface condition, the largest 
single contamination patch centered at 700 feet past the exit entrance (full exit width) is 
the following: (see Maximum Ground Coefficient Friction for MD-1 1 Diagram, figure 6.4) 

50 foot long patch of Textured Concrete Flooded 

OR 

15 foot long patch of Ice (0.02*Dry) 

Friction Patches affect Exit Prediction Logic 


The exit prediction logic using speed dependent measured friction along the runway at 250 
foot intervals appears to be quite reliable. However, studying varying patches of friction has 
shown that the prediction can be false if the measured runway friction has errors due to 
friction patches not included in the measurements. Either the friction interval should be as 
small as the smallest reasonable patch size (5-15 feet) or the measured friction should be an 
average of the friction over the runway interval (250 feet) it represents. The current exit 
prediction logic makes a few more incorrect predictions as available friction becomes more 
critical with lower friction conditions. 
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LITERATURE SEARCH 


Introduction 


The ROTO system was designed to provide safe runway operation on dry and wet runway 
surfaces. The problems of operation for conditions or runways being contaminated with 
water (flooded), ice, snow-packed, or slush have been identified and are discussed in this 
section. 

A goal of the research described herein was to develop a ROTO system to decelerate and 
steer the aircraft onto high-speed runway exits and onto taxiways at entrance speeds to the 
high-speed exit up to 70 knots. A suitable guidance and control (G&C) law to automatically 
steer a transport aircraft along runway centerline, onto a high-speed runway exit, and finally 
along a runway taxiway was to be developed. The system was also to be capable of 
decelerating the aircraft by means of autobraking and auto-reverse thrust to the desired turn- 
off and taxi speeds. 

Adequate tire-runway interface friction for braking must be available for the system to effect 
necessary longitudinal deceleration and lateral control. This section addresses the tire-runway 
interface friction. 

Operations 

General 

The runway friction during the high-speed turnoff is critical; therefore, it is 
recommended that the runway in the vicinity of the turn be grooved. Further, the 
grooving should be in a “diamond” pattern. This type pattern will improve the friction 
coefficient in all directions. 

During the course of this study, CUSHION CUT, a manufacturer of pavement grinding 
and grooving equipment, in near-by Torrance, CA, was visited. It was observed that all 
passage-ways in their factory area were grooved. At the intersections there was cross- 
grooving. From walking on the floor in shoes with hard leather soles it was noticed, 
where the grooving was longitudinal, the lateral friction was high; whereas, in passage 
way intersections, with a “diamond” groove pattern, the friction was high in all 
directions. 

This small example is believed to illustrate that cross grooving in the vicinity of the turn 
would resist the tendency of the aircraft to skid laterally when the centrifugal forces are 
highest. 
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If moisture is on the pavement, grooving will also prevent reverted rubber hydroplaning 
which may be started prior to the aircraft’s reaching the point of turning. In this type 
of hydroplaning a pad of steam develops under the tire eliminating virtually all friction 
between the main tires and the pavement preventing braking and directional control; 
therefore it must be avoided. 

Wet and Dry Runway 

In beginning this discussion it is worthwhile to classify the “wet” runway as one which 
is well soaked but with no measurable amounts of water. 

An international group consisting of the Federal Aviation Administration (FAA), Joint 
Aviation authorities (JAA) (European equivalent of FAA), Civil Aviation Authority 
(CAA) (British equivalent of FAA), Transport Canada (TC), National Aeronautics and 
Space Administration (NASA), and Manufacturers: Douglas, Boeing, Airbus, 
Gulfstream, Alenia, & Aerospatiale have been endeavoring to develop civil regulations 
for take off field lengths that account for degraded stopping performance on wet 
runways. To derive a conservative wet to dry friction ratio, data from some B737 tests 
demonstrated wet runway landing certification results, in which the average airplane 
braking coefficient on wet, grooved pavement exceeded 90% of the demonstrated dry 
runway value. This is illustrated by figure 6. 1 . Additionally, selected pages from 
NASA report SP-5073 (Reference 3 is shown by figures 6.2 and 6.3). These illustrate 
various runway surfaces and degrees of wetness and flooding. On each of these is a 
curve of 70% of dry friction coefficient. It may be seen that this curve shows adequate 
friction margin for all depths of wetness when the surfaces are grooved. 

From this discussion, the conclusion may be drawn that for the WET and DRY 
RUNWAY, when grooved and adequate coefficient of friction margin exists, no problem 
will exist in accomplishing the high-speed rollout and turnoff at 70 knots. 

Contaminated Runway and High-Speed Exit 

Advisory Circular No. 91-6B (Ref. 4) defines a contaminated runway as one with more 
than 25% of the required field length, within the width being used, covered by standing 
water or slush more than 0. 125 inch (3.2 mm) deep, or that has accumulation of snow or 
ice. Other situations would be appropriate to say it is contaminated if there is an 
accumulation in critical areas such as the high speed for lift-off, or the high-speed exit. 
The runway friction is as an average value. An average value is suitable for the runway, 
but not for the exit. The specific exit friction must be determined. For, if the aircraft 
encounters a slick area while turning, it may skid out of control. 
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Normal Deceleration 


An assumption for ROTO operations is that the runway will have a first, second, and 
third high-speed exit. Nominally on dry or wet runways, the airplane would decelerate 
sufficiently to approach the first or second exit at 70 knots or less. The highspeed exit 
would then be accomplished. 

Low Deceleration 

If the airplane could not decelerate to 70 knots at the first exit because of contamination 
and a low friction level, the first exit would be aborted. If the deceleration is low enough 
that 70 knots is exceeded at the second and third exits, these would be aborted also. 

The aircraft would have to continue to the end of the runway before turnoff at a low 
speed. 

If an exit is aborted, a closely-spaced following aircraft would likely have to execute a 
go-around, and the ROTO sequence would need to be modified. 

Run way/Taxi way Surface Condition Sensors 

It is recommended that Runway Surface Condition Sensors in accordance with (Ref. (5) 
be installed in the runway and highspeed turnoff at locations and in quantities as 
necessary to monitor critical areas. The highspeed turnoffs would likely be the most 
critical locations requiring the greatest number of sensors. Otherwise, the general locale 
of the airport - its proximity to mountains, bodies of water, plains, etc.-, would dictate 
their need and location on the runway. 

Signals from the surface condition sensors, embedded in the exit surface, would be 
monitored by the airport control tower. Evaluation of the sensors measurements would 
allow the ROTO system to determine whether or not an automatic high-speed turnoff 
could be accomplished. 

Processing of the surface sensor measurements could indicate the acceptable speeds for 
use of the high-speed exit. As may be seen on friction curves in figure 6.4, at this speed 
adequate friction is available except for snow conditions to prevent lateral skidding of 
the airplane. 

Figure 6.4 also shows ice to have a coefficient of approximately 0.05; therefore, 
movement of the airplane at the end of the runway would likely be accomplished with 
an airplane tow-bar and tug. 
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In addition to data from the surface sensors, friction information about the runways and 
exits should be determined as necessary with continuous friction measuring equipment 
as described in reference 1 . 

Discussion of Contamination 


General 

Under this heading the problems for potential operation on runways which are 
contaminated by flooding, ice, packed snow, or slush are discussed. These are the basic 
natural contaminants, which may be present in various forms, primarily due to 
temperature. The biggest problem arises from any one of these contaminates if the 
quantity is so great that it cannot be controlled. This results in closure of an airport, 
and if the airport is a major hub, airports throughout the nation will be affected. Early 
in January 1996, this was the case with many airports in the Eastern and Northeastern 
part of the United States. This was extreme, but lesser extremes can and will affect the 
operation of a High-Speed Rollout and Turnoff Guidance and Control System. The 
effects of these extremes can be reduced and the system can be back in operation sooner 
when certain precautions are made. Grooving the runway and exits is a highly 
important consideration. 

The basic effect of these extreme conditions is the closure of the airport or a particular 
runway . The decision for closure can likely be made with more certainty than whether 
or not an exit is capable of supporting a high-speed turnoff Several factors affect the 
precise determination of the friction level -- the particular contaminate; its roughness 
and hardness, the condition of the runway surface, macro and/or micro surface texture; 
temperature, condition of the airplane tires; calibration of the friction measuring devices; 
their correlation with the aircraft being considered; and more. The FFA (The 
Aeronautical Research Institute of Sweden) has made extensive studies concerning 
snow, ice and slush on runways (Reference 6). One result showed that runway 
temperature around 0 degrees C in combination with precipitation very often caused 
great discrepancies between experienced and measured runway friction. 

An average friction value for the highspeed turnoff is inadequate. The minimum 
available friction value for a high-speed exit must be known. If a higher value is required 
at any point during the turn for the desired turnoff speed, loss of lateral control of the 
airplane could result. Thus, an average friction value for the high-speed turnoff is 
inadequate. 

Flooding 

Reference 7 summarized the research done relative to flooding during rainstorms, 
hydroplaning, identification of slippery runways, progress in developing 
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antihydroplaning procedures, and rubber deposit removal. Grooving is extolled - 
particularly grooves made with saws having diamond chips embedded in their cutting 
edge. 

The exit speed of 70 knots is below the full hydroplaning speed of the MD-1 1; 
however, at this speed there is a wedge of water under a portion of the tire. If the 
pavement is grooved, the area of the wedge is reduced and more of the tire is in contact 
with the pavement even though the rate of rain fall is sufficient to cause complete 
flooding. Thus, some data supports the capability of making a high-speed exit on a 
grooved flooded runway. This capability may depend on the details of all the variables 
involved. Notwithstanding, based on the figure 6.4 friction curves the ROTO 
Operation at 70 knots cannot be accomplished on flooded textured concrete. 

Pack-Snow, Ice and Slush 

Advisory Circular No. 91-6B gives guidelines, and recommendations concerning takeoff 
and landing performance. It indicates that provisions should be for test measurements 
to establish the runway braking coefficients for standing water, packed snow, slush and 
ice. If such data is not available, appropriate data supplied by the manufacturer should 
be used. For the aircraft being considered in this study, test data for all conditions are 
not available, but training simulators data has been derived figure 6.4. This figure 
indicates the coefficient for packed-snow is marginal, for slush, the friction is high; but 
for ice, the coefficient is too low to consider. 

These data correlate reasonably well with data from Ref. (8). This reference reports the 
Aircraft and Ground Vehicle Friction Correlation Test Results Obtained Under Winter 
Runway Conditions During Joint FAA/NASA Runway Friction Program. The aircraft 
used for these tests were the B-727 and the B-737. The mu effective values recorded 
for packed snow at 70 knots were between 0.15 and 0.21. The values for solid ice were 
0.03 to 0.05 at 70 knots. There was no slush tested during this program. 

Procedures for Minimizing the Affects of Contaminates 

Weather in various parts of the world can be highly variable. Rain squalls can cause one 
part of a runway to be flooded for some time while the other end is completely dry. 
Freezing rain will stick to the pavement if it is already at or below freezing 
temperatures. High gusting wind can make precise control of the aircraft difficult. The 
snow fall can be so great that it exceeds the capability of the maintenance equipment to 
remove it adequately. Aircraft cannot land or take off from the airport. In extreme 
circumstances it must be closed. If precipitation rates are not too high and the 
temperature is not too low, then with judicious maintenance the airport surfaces can be 
kept free of contaminate such that they are not worse then wet conditions. 
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Two important advisory circulars are AC No. 150/5320-12B (Ref. 9), and AC No. 
150/5200-30 (Ref. 10). The first gives guidelines and procedures for design and 
construction of skid-resistant pavements; pavement evaluation, and maintenance of 
high-skid-resistant pavements. This circular covers both the Asphaltic Concrete and 
Portland Cement Concrete. There is a good section on grooving for both types of 
concrete. Information is given relative to friction deterioration, evaluation and 
frequency of survey. The need for continuous friction measuring equipment (CFME) 
and general specifications are given. 

The second advisory circular provides information to assist the airport owner/and 
operators in the development of acceptable airport snow and ice control programs. 
Chapters are given on personnel and equipment organization for Winter Operations on 
Airports; Airport Snow Removal Equipment — Equipment selection, Services, and 
minimum equipment. Snow and Ice Removal Procedures — Snow control procedures, 
Mechanical Methods and Chemicals for control and removal of snow and ice, and 
chemical affects on friction; also, use of abrasives and Runway Friction Improvements. 

Neither of these circulars are strictly directed toward, or provide specific instructions 
relative to the High-speed Turnoff, but they serve to show the great amount of 
organization, equipment and care necessary to assure winter safety. The requirements 
for ROTO are more critical than addressed by the circulars and will require more study 
to assure that all facets are recognized and satisfactory solutions are developed. 

Summary and Conclusions 

Problems which will be encountered during the ROTO Operation with contamination of 
flooded water, ice snow-packed, and slush are discussed. Grooving of the pavement surface 
will enhance the friction, particularly for contaminated conditions. A diamond groove pattern 
will tend to cause improvement in all directions and is recommended in the area of the exit 
turnoff. 


Dry and wet (well soaked, but with no measurable water) runway/taxiway appear to be 
satisfactory for accomplishing the high-speed rollout and turnoff at 70 knots. 

Possible scenarios are discussed with varying amounts of contamination and location on the 
runway & exits and the use of surface condition sensors. These devices strategically placed 
can be monitored to indicate whether or not a safe turnoff can or cannot be made at a 
particular exit or subsequent exit locations. 

The runway friction coefficient must be measured at regular intervals along its length. This 
data would be used by approaching aircraft. 
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Precipitation rates cannot be controlled. If the airports’ available maintenance/cleaning 
equipment can cope with the rate, the pavement will be no worse than wet for most 
contaminates except for deluges of rainfall, snowfall and freezing rain. 

Highspeed exits may be made for dry and wet conditions, slush, and for “flooded” if the 
pavement surfaced is grooved and depth of water is not too great. 

Highspeed exits should not be made on ice or packed-snow. The latter appears to be marginal 
and should not be considered until further research and/or procedures are developed. 

Recommendations 


The goal of the ROTO concept is to provide a capability for landing airplanes every 50 seconds. 
Undoubtedly there will be some weather and contamination conditions which will prevent this desire 
from being met. The following considerations may help in meeting this criteria and/or clean-up of the 
runway and turnoff during or after severe weather conditions. 

On various runway surfaces the depth of tire tread may make a considerable difference in the tire- 
runway interface friction. A method of determining and reporting the tire tread conditions would be 
beneficial. 

Reverted rubber hydroplaning is hazardous. Precise reasons for this occurrence should be determined 
and measures taken to prevent. 

When the 70 knot turn off cannot be met, a procedure for turnoff at lower speeds should be 
determined. This will reduce the runway occupancy time to less than the airplane being taxied to the 
very end of the runway. 

There is a potential for the intersection of the runway and the highspeed turnoff, when the turnoff 
has a large radius, to have a large flat area which would drain poorly. Cross sections of such 
intersections should be studied to obtain optimum drainage. 

Timing for the cleaning of the runway and taxiway after the subsidence of severe weather 
may be critical, so, the coordination for removal of contaminant and the maintenance of the 
airport/runway must be precise. 
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7.0 ROTO DESIGN 


All of the ROTO control laws have update rates at 20 Hz. In this study, all control 
calculations that create the rate of an input are suspended if the input data is not fresh or is 
dropped out. The navigation inputs are updated at 10 Hz. All control laws appeared to be 
equally suitable for both MD- 1 1 and MD-8 1 type aircraft. 

AUTOBRAKING CONTROL LAW 


The autobraking control law provides a variable deceleration command to the autobrake 
system. The purpose is to decelerate the aircraft to the selected ROTO exit’s entrance speed 
and then continue to decelerate on the ROTO exit so that the aircraft may stop prior to the 
taxiway if desired. The deceleration command is limited to 9 ft/sec 2 (heavy braking) and rate 
limited so as to minimize longitudinal jerk. Some MD-1 1 simulation runs approached 9 
ft/sec 2 during a portion of the runway. Inputs to the autobraking control law are the aircraft 
ground speed, aircraft X,Y position, selected ROTO exit entrance X position, selected ROTO 
exit entrance ground speed (70 knots used in this study), taxiway entrance X,Y position and 
taxiway entrance ground speed. The top level diagram of the autobraking control law is 
illustrated in Appendix A on page A-88, with lower diagrams thereafter. A simplified 
representation is shown in figure 7.9. 

This study assumed that the pilot engages reverse thrust and speed brakes as soon as 
possible after touchdown, which will minimize the required braking. Auto-Reverse Thrust is 
discussed in the next section. The autobrake system is referred to in reference 1 or in 
Appendix A on page 125. A commanded deceleration results in brake pressure. The 
actuation system is identical to the existing MD-1 1 autobrake system except that instead of 3 
available settings (high, medium, low), the actuation model allows the setting to be 
continuously varied by the autobraking control law. As in the existing system, the model 
delays all braking until after spoiler deployment to ensure adequate vertical loading of the 
main gear. Full braking ramp rates and gains are only available after nose wheel touchdown, 
which is approximately 6 seconds after main gear touchdown. 

In order to minimize runway occupancy time (ROT), it is desired to delay the onset of 
braking until a ROT less than 53 seconds is assured. Safety issues might argue against this 
practice. ROT begins counting when the aircraft is at the runway threshold. The logic to 
begin braking is as follows: 

1 . Estimated final ROT is less than 53 seconds. OR 

2. Estimated final ROT is increasing, ignore reverse thrust spool up time period. OR 

3. Estimated final deceleration is greater than 6-8 ft/sec 2 (logic dependent). OR 

4. A 53 second ROT requires a constant deceleration greater than 7.5 ft/sec 2 . 
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Logic item 1 minimizes the deceleration magnitude once it is guaranteed that the estimated 
ROT is under our desired value. For logic item 2, if estimated ROT is increasing we assume 
that it is best to begin braking for the current exit rather than consider any later exit, even 
though the other logic may show that the estimated ROT is greater than our desired value. 
This logic should use a steady state estimated ROT, so we ignore the time during reverse 
thrust spool up. This study assumed reverse thrust spool up has ended 6 seconds after main 
gear touchdown. Logic item 3 uses current deceleration and deceleration rate to estimate a 
final deceleration. Logic item 4 attempts to keep a constant required deceleration under 7.5 
ft/sec 2 for the desired ROT value. 

When braking begins, a linearly decreasing speed profile versus the distance to the selected 
exit is created. The deceleration command is created by a PI controller, whose input is the 
aircraft’s ground speed minus the desired speed at that location on the runway. This speed 
versus runway location profile is more easily modeled by the exit prediction logic described in 
a later section, rather than an initial deceleration method (reference 1 ) which decelerated the 
aircraft to reach a final speed at the exit entrance only. Because there is a target speed all 
along the runway, intermittent runway low friction patches are more quickly corrected for. 
The speed profile has a distance buffer (~100 feet) just prior to the selected exit. 

The selected exit is aborted, up until 250 feet prior to the exit, if the aircraft’s current required 
deceleration or estimated final deceleration is greater than 9 ft/sec 2 . At the exit entrance the 
exit would be aborted if the aircraft’s ground speed is greater than 2 knots above the exit 
entrance speed or that the aircraft has a lateral displacement from the runway centerline of X 
feet (needs to be investigated). MD- 1 1 simulations on a wet runway with a steady 1 5 knot 
crosswind did not see a lateral deviation too great to cause an exit steering problem. If an exit 
is aborted, the aircraft would revert to non-ROTO. More discussion is required to decide 
whether pilots would then again want to again select auto ROTO for the next available 
ROTO exit. 

The exit prediction logic, prior to touchdown, attempts to minimize the occurrence of 
aborting an exit on the runway by the deceleration exit-abort logic in the preceding paragraph. 
Before the exit prediction logic was implemented, the affect of the preceding paragraph’s exit 
abort logic was: 

1 • The aircraft s current required deceleration would either abort an exit(s) immediately at 
touchdown or during the last third of the distance to a selected exit. 

2. The aircraft’s estimated final deceleration generally would abort an exit during the first 
third of the distance to a selected exit. 

ROT is minimized by aborting appropriate exits as soon as possible, preferably prior to 
touchdown. 
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AUTO-REVERSE THRUST CONTROL LAW 


Reverse thrust is needed for operations in low friction runway conditions and to minimize 
brake usage. An assumption for this study was that the pilot moves the throttle levers 
through the pedestal inter-locks to engage reverse thrust soon after touchdown and stows 
reverse thrust (or at a minimum sets it to idle) at 70 knots ground speed. Currently pilots 
generally stow reverse thrust based on airspeed. It is suggested that pilots stow reverse 
thrust just prior to entering the exit to negate that task while on the exit. 

When reverse thrust is engaged, the auto-reverse thrust control law commands the auto- 
throttle servo to drive the throttle levers to achieve the desired reverse thrust. The control 
law adjusts reverse thrust so that the symmetric brake pressure stops between 20 - 30 % of 
maximum supply pressure during the ROTO deceleration phase. Because the engines have 
slower dynamics than the brakes, the brakes provide the longitudinal damping. The reverse 
thrust is allowed to vary between idle and maximum reverse thrust (Note: These negative 
thrust limits become more positive as airspeed decreases). When the deceleration command is 
zero, symmetric brake pressure decreases to return pressure. 

The auto-reverse thrust command is integrated to slowly decrease reverse thrust to idle if 
brake pressure is below 20 % or the aircraft’s ground speed is below the desired speed 
profile. There is a command deadband when the brake pressure is between 20 - 30 %. The 
command is integrated even more slowly to increase reverse thrust to maximum if brake 
pressure is above 30 %. The figure in Appendix A on page A-87 illustrates the auto-reverse 
thrust control law. 
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STEERING CONTROL LAW 


The PID steering control law used for this study started with a basic auto rollout control law 
similar to localizer type control laws. Lateral centerline deviation is the main command error 
with lateral centerline deviation rate and yaw rate as damping feedbacks. On the ROTO exit, 
exit centerline radius and radius rate were also added as feedbacks to lead the aircraft into the 
curve. Aircraft ground speed was used for gain scheduling. The output of the control law is 
the position limited rudder command. . The top level diagram of the steering control law is 
illustrated in Appendix A on page A- 106, with lower diagrams thereafter. A simplified 
representation is shown in figure 7. 10. The nose gear is geared to the rudder through the 
rudder pedals and is limited to +/- 8 degrees. Navigation inputs (both the aircraft position 
and the data stored for path centerline) were converted to a X,Y reference frame for lateral 
deviation calculations. Nose gear loading was maximized by co mmandin g a positive down 
elevator, scheduled with ground speed and aircraft CG. 

For typical exit path geometry, lateral jerk normally occurs as the aircraft begins to follow the 
centerline of the exit because the physics of the path geometry require a step change in lateral 
acceleration to follow the exit path centerline. Gradual spiral exit (reference 2) entrances help 
to minimize this effect. If the exit has a constant radius entrance, the control law can contain 
a feature to smoothly increase the lateral acceleration to nominal lateral acceleration as defined 
by ROTO exit radius and desired exit speed. The lateral acceleration is increased at maximum 
allowable lateral jerk. This procedure is begun at a distance prior to the exit based on the 
aircraft ground speed and how many seconds are required to ramp up the lateral acceleration. 
During the pre-exit lateral acceleration, the lateral deviation of the aircraft changes less than a 
foot due to the commanded lateral acceleration. 

With further steering control law development, auto asymmetric braking, geared to the rudder 
command to augment steenng, was found in this study to be unnecessary. This study 
showed that nose gear steering with no more than 2 degrees of hysteresis is suitable for an aft 
CG MD-1 1 on a wet runway. It is suggested that the tiller be used to compensate the nose 
wheel such that the MD-1 1 auto ROTO nose gear hysteresis will be effectively reduced to 2 
degrees. This suggestion for nose gear steering was implemented for this study in a tiller 
control law. The tiller control law subtracts measured rudder pedal angle from the desired 
nose gear (rudder) command. After washing out this signal, it is used to command a “cable 
grabber” servo which is connected to the tiller-to-nose gear cable. The tiller handle is back 
driven and the servo is commanding the nose wheel actuator through the tiller cables. This 
signal is mechanically summed with the rudder pedal nose wheel command. The tiller control 
law is shown in the figure in Appendix A on page A-l 13. It is believed that the tiller control 
law would be less costly than a fly-by- wire nose gear to decrease hysteresis of the baseline 
MD-1 1 auto nose wheel steering. 

In auto ROTO the tiller control law activates prior to turning onto the exit entrance. The 
tiller control law is not active for manual ROTO. In manual ROTO, the pilot follows the 
flight director guidance with the yaw damper on. 
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AUTO ASYMMETRIC BRAKING CONTROL LAW 


Reference 1 refers to the need for asymmetric braking to assist in steering an aft C.G. MD-1 1 
on a wet high speed ROTO exit. Development of the baseline steering control law has now 
alleviated the need for this feature, as well as a 1 -degree hysteresis, fly-by-wire nose gear 
(assuming auto ROTO tiller control law requirement). However, asymmetric braking would 
improve steering if it were available. All simulation runs described in this report do not 
include auto asymmetric braking, nor a fly-by-wire nose gear. The developed auto 
asymmetric braking control law is illustrated in the figure in Appendix A on page A-94. 

The auto asymmetric brake command is geared to the rudder command and summed with the 
symmetric braking command at the brake. The asymmetric command passes through a 
deadband to minimize deceleration braking. An opposite command is sent to the opposing 
brake. If the command is in the direction to drive a brake’s pressure below its return 
pressure, the command is doubled to the opposing brake. 

The command is rate limited to minimize longitudinal jerk. Longitudinal jerk due to the brake 
pressure threshold cannot be helped, except to not let the brake pressure fall below the 
threshold during the asymmetric braking phase. Whenever the command is zero, a temporary 
command is created to drive the brake pressure asymmetry back to zero. Because the brakes 
act as integrators, if this were not done, some type of asymmetry would persist thereafter 
and fight the commanded rudder and nose gear. 
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EXIT PREDICTION LOGIC 


In order to minimize runway occupancy time by controlled deceleration, it is desired to 
predict which available ROTO exits the aircraft is capable of using at touchdown. For ROTO 
pilot arming/work load reasons, this prediction would most likely occur up to a half minute 
prior to touchdown if touchdown parameter predictions are accurate enough to allow this. 
Whether a ROTO exit is selected by the pilot or by an exit prediction algorithm; selecting too 
early an exit would cause the exit to be aborted, causing the aircraft to coast to the next exit. 
Selecting a late exit would increase runway occupancy time above what necessary. Both of 
these occurrences may cause the following aircraft to go-around. 

The exit prediction logic has the following predicted/estimated inputs: touchdown location, 
touchdown ground speed, aircraft weight, aircraft CG, aircraft drag characteristics and aircraft 
thrust versus airspeed/time profiles. Estimated steady airport winds and speed dependent 
runway friction along its length are also required by the algorithm. It may be possible to 
predict autoland touchdown location more easily than ILS manual or non-ILS landings due to 
the flare segment. MD- 1 1 ROTO simulations have recommended the following MD- 1 1 
predicted value accuracies. The Exclusive Accuracy column is the maximum allowed error for 
a predicted value with that error being tested by itself. The Combined Accuracy column is 
the maximum allowed error for a predicted value with all of the errors being tested together. 


Predicted Value 


Exclus ive Accuracy Relative to Actual Combined Accuracy 


Touchdown location 
Touchdown ground speed 
Touchdown weight 
Runway friction coefficient 


0 to +300 feet 
0 to +3 knots 
0 to +100,000 lb. 
-0.03 to 0 


0 to +100 feet 
0 to 1.0 knots 
0 to 30,000 lb. 
-0.01 to 0 


It appears that the current modeled MD-1 1 dispersion and exit locations only allow the 
prediction errors to lie on the conservative side. If the predicted touchdown location error 
were -100 feet, the exit prediction logic would perhaps accept too early of an exit which 
might be eventually aborted. If the predicted touchdown location error were +400 feet, the 
exit prediction logic would perhaps accept to late of an exit than necessary which would 
unnecessarily increase the runway occupancy time. 

These prediction accuracies are dependent on the simulation accuracy, the number and 
location of ROTO exits and aircraft type. It may not yet be technically feasible to achieve 
these accuracies necessary to correctly predict ROTO exits, which is essential to minimize 
following aircraft go-arounds which are currently at about 1 per 1000 landings. The 
prediction errors would decrease as the prediction is made nearer to touchdown. 
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These MD-1 1 simulations assumed speed dependent friction was measured along the runway 
approximately 250 feet apart. It was found that this friction measurement should 
represent/average the friction value of the entire 250 feet. If there are friction patches on the 
runway and the friction measurement is in or out of a patch, these non-representative 
measurements may accumulate to cause a false exit prediction. 

The exit prediction algorithm attempts to simulate the aircraft’s onset of deceleration and 
deceleration (due to brakes, thrust, drag, speed brakes) down the runway until the aircraft CG 
reaches the target exit entrance. The logic to begin aircraft braking is referenced in the auto- 
braking section above. The exit prediction simulation steps down the runway in terms of 
distance rather than time. Because the actual aircraft attempts to decelerate according to a 
linear speed profile versus distance rather than target a final speed at the exit entrance, it is 
more easily modeled by the exit prediction algorithm. 

When the simulated aircraft reaches the predicted exit, if the required friction is greater than 
the available friction or a deceleration greater than 9 ft/sec 2 is required, the predicted exit is 
aborted. The prediction algorithm then repeats until a suitable available exit is found. 

The exit prediction logic path is illustrated in figure 7. 1 and MATLAB script code in 
Appendix A on page A-27. 
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FLIGHT DIRECTOR MANUAL ROTO GUIDANCE 


Manual ROTO steering and deceleration flight director guidance control laws were not 
designed as part of this study. The flight director commands would most likely be a filtered 
output of the auto ROTO steering and deceleration control law’s forward paths. The yaw 
damper should be active with manual ROTO steering. The HUD display provides both 
command and raw data guidance. 
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CURRENT ROTO OPERATIONS 


Some U.S. airports, e.g. Dallas-Ft. Worth, under flight crew discretion conduct manual high- 
speed ROTO under daylight VMC conditions with no surface contamination on the exits. 
This applies to both narrow and wide body aircraft on 30 degree exits at exit entrance ground 
speeds up to 60 mph. 

The current runway clearance definition which would apply to continuous ROTO operation 
is as follows. The runway would be clear for a following aircraft at the runway threshold, if 
all preceding aircraft have cleared the runway and are past their exit hold lines (figure 7. 1 1 ) or 
are in the process of rolling past their exit hold lines without obstruction. The runway would 
not be clear if a preceding aircraft has stopped on an exit prior to an exit hold line. 

If an aircraft stops on an exit past the hold lines, the taxiway is open, but that exit is not 
available as a ROTO exit for following aircraft. The runway would not be clear for a 
following aircraft that required that occupied ROTO exit, unless the following aircraft 
performs a NON-ROTO landing. 
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ROTO MODE LOGIC 


Figure 7.2 illustrates the ROTO modes (states) and the transitions between those modes. 
For any mode, the arrow originating from a darkened circle is the default transition. In 
approach, ROTO begin in the NON ROTO state. After 

ATC has provided the expect ROTO clearance AND 
ROTO exit prediction software has identified a ROTO exit 

, ROTO transitions into the ROTO READY mode. The pilot may arm ROTO and an exit 
causing the transition into the AUTO ROTO SELECTED mode by pressing the ROTO 
button on the glareshield. AUTO ROTO provides automatic steering and deceleration on the 
runway and onto and through the ROTO exit, while the pilot monitors out-the- window 
through the HUD and other displays. 

The pilot may transition between the READY, AUTO ROTO and MANUAL ROTO 
(manual steering with guidance/while maintaining auto deceleration) modes by repeatedly 
pressing the tri-state ROTO button on the glareshield. Additionally, after touchdown the 
pilot may also press the yoke autopilot disengage button to select MANUAL ROTO. 

After touchdown in AUTO ROTO SELECTED mode, if the pilot overrides the rudder 
pedals, ROTO would transition to MANUAL ROTO (manual steering with guidance/auto 
deceleration) mode. After touchdown in AUTO or MANUAL ROTO, if the pilot overrides 
the brake pedals or throttle levers, ROTO would transition to FULL MANUAL ROTO 
(manual steering with guidance/manual deceleration with guidance). 

At any time, AUTO or MANUAL ROTO would revert to NON ROTO if 
Pilot deselects ROTO OR 
ATC rescinds ROTO Clearance OR 
Aircraft is in GO- AROUND OR 
ROTO exit is aborted OR 

Any mechanical/software ROTO function has failed on the runway prior to the exit. 
The failure must be positively indicated. 

Prior to touchdown the pilot may select NON ROTO by pressing the yoke autopilot 
disconnect button. After touchdown the pilot may toggle the glareshield ROTO button to 
the READY state. 

ROTO would remain in its current mode if any mechanical/software ROTO function fails on 
the exit. A failure annunciation would notify the pilot of the condition, whereupon he/she is 
free to override the controls and manually control using ROTO guidance 

Under certain timing conditions, if an exit is aborted (ROTO mode reverts to NON-ROTO), 
there may be enough remaining runway occupancy time to transition again from NON-ROTO 
to AUTO ROTO by selecting a second available ROTO exit, then exiting that second exit 
while still keeping the runway occupancy time under 50 seconds, thereby not requiring a 
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close following aircraft to GO-AROUND. This topic needs further discussion and pilot 
input. 
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ROTO OPERATIONS AND CREW INTERFACE 


The ROTO crew interface concept was developed in response to anticipated flight crew and 
Air Traffic Control (ATC) operational requirements associated with ROTO's implementation 
in terminal area operations. More specifically, development of the crew interface concept for 
ROTO operation was accomplished in two phases. In conjunction with a McDonnell 
Douglas engineering test pilot and two commercial airline pilots, a targeted operational 
analysis was conducted articulating approach, landing, and roll out functions required of the 
crew and ATC. Based on this analysis, a design effort was performed — again with 
substantial pilot contribution and review - yielding the specification of control and display 
system elements needed to perform these required activities. 


Operational Analysis 


To conduct a more substantive operational analysis, crew system design personnel and the 
engineering test pilot familiarized themselves with the major functional capabilities and 
system logic parameters of the ROTO system: Relationships between ROTO-specific 
system calculations and performance, and potential modifications to normal guidance, brake, 
and reverse thrust control; automatic and manual modes; and deceleration and steering 
profiles. Of specific interest were estimates of ROTO calculation times, crew-ROTO 
coordination in the application of thrust-reversing and braking, and the requirement for out- 
the- window monitoring of ROTO progress. 

After acquiring general familiarity with ROTO system performance, a representative 
functional timeline was developed to delineate how ROTO operation should be integrated 
with standard terminal area arrival, approach, and landing procedures. This timeline, 
presented in Appendix A on page A- 139, identifies the necessary communication events 
between the crew and ATC, and indicates the probable sequencing and timing of crew and 
ROTO system activities. The timeline was particularly useful for articulating the relatively 
critical timing constraints associated with ATC's clearance for ROTO and the subsequent 
execution of ROTO-directed landing, roll out, and turn off. (It should be noted that, in the 
timeline, ATC maintains ~50 second separation between ROTO-controlled landings by 
issuing the ROTO clearance for the aircraft under discussion immediately after a preceding 
aircraft has cleared the exitway needed for the landing aircraft - in the last five to ten seconds 
of flight before touch down. This proposed procedure must be carefully evaluated in 
subsequent ROTO operational analyses because of concerns about workload, crew situation 
awareness, and the actual feasibility of ATC issuing such clearances in such small and critical 
'windows' of time. Alternative procedures, such as earlier clearance issuance with possible 
(last minute) recision, would need to be explored in on-line piloted simulation activities.) 
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ROTO Control and Display Elements 


From this operational analysis, control and display elements were designed and their 
operations defined. Following a design philosophy of minimal modification of existing 
procedures, and correspondingly minimal modifications of control and display elements, the 
ROTO system as shown in Figure 7.3 comprises several components: 

1 . A modified Head Up Display. 

2. A glareshield control head (including ARM and READY settings) and annunciation 
panel. 

3. A Primary Flight Display containing a backup (head-down) representation of ROTO- 
related information displayed on the HUD. 

4. A Navigation Display (ND) showing a plan view representation of the aircraft's 
progress on the ROTO deceleration profile and steering track. 

5. A slightly modified Flight Management System page for specifying and executing 
ROTO calculations/estimates (to prepare ROTO for initiation after clearance). 

6. An optional data link system for receiving ROTO-relevant airport information, and for 
crew acknowledgment of ROTO's clearance. 

7. A modified Autopilot Disengage switch on the yoke for ground-operational de- 
selection of Automatic ROTO steering (thus releasing steering control to the crew for 
Manual-mode ROTO operation). 

8. A "ROTO" setting on the Autobrake Control panel. 


Owing to the limited and relatively preliminary scope of this activity, the design effort was 
focused exclusively on the HUD and glareshield components -- elements critical to ROTO's 
implementation, and requiring relatively extensive operational/procedural analysis. 
Modifications to the FMS page(s), PFD and ND formats, data link, and Autopilot Disengage 
switch were viewed as straightforward, and were, therefore, not developed in this design 
effort. 

As can be seen in Figure 7.4 the HUD contains a number of ROTO operation-specific 
display elements. Positioned above the horizon line and in the center of the HUD field are 
the ROTO and exit number designation and status elements. The word "ROTO" is followed 
by the system-selected exit (in this case, " 1 "), and both are enframed by a box symbol. The 
box is written in a solid line and is labeled, in this example, with the designation "AUTO," 
indicating ROTO's present status (i.e., being armed to execute a ROTO procedure in which 
steering is automatically controlled by ROTO). Before the crew had armed ROTO, the box 
was shown in a dashed line, and the system’s "READY" status was displayed. This 
information is replicated on the glareshield (described subsequently). 
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To the left of the ROTO/ exit number status indication elements, and also above the HUD 
horizon line, is the Airport and Active Runway designation. This information (also repeated 
on the glare-shield) is simply provided as a head-up verification that the FMS-specified 
selection is indeed correct (this feature is a useful safety enhancement element, especially in 
cases of reduced visibility, problematic ATC communications, areas with closely co-located 
airports, and for airports with multiple runways). 

Two symbols commonly employed (in various forms) on ‘glass cockpit’ PFDs and HUDs « 
the Deceleration Cue (a ">" symbol positioned to the left of the Flight Path Vector symbol) 
and the Speed Deviation Indicator (a deviation bar rising from, or descending from, the left 
wing of the Flight Path Vector symbol) - have been co-opted to support guidance and 
control capabilities associated with following the ROTO system-calculated speed schedule. 
As in standard usage, the Deceleration Cue drifts below alignment with the Flight Path Vector 
symbol's left wing, commensurate with the aircraft's current (ROTO-commanded) 
deceleration. Precision guidance and/or control is provided by the Speed Deviation Indicator 
which, in turn, is driven by the autobraking control law. Length of the deviation bar indicates 
magnitude of the instantaneous speed deviation (from directed speed schedule), and direction 
(up or down) of the bar indicates instantaneous speed error, with an upward extended bar 
indicating that the instantaneous speed is above commanded and a downward extended bar 
indicating below commanded. 

The Flight Path Vector symbol, a circle with two down angled "wings" (one to the circle's left 
and one to its right), operates in a standard HUD-operational manner, indicating the 
instantaneous flight path of the aircraft. ROTO operation extends the use of this symbol to 
the indication of the ground path vector. This information facilitates runway and exitway 
centerline tracking, and assists in exit- and taxiway turn guidance (especially in reduced 
visibility). The Flight Director symbol (a circle sized to be clearly distinguishable inside the 
Flight/Ground Path Vector circle) is also used as guidance for pursuit of the desired ROTO 
path centerline, and is of critical importance for precision tracking in the complex ground turn 
profiles (spiral profiles) prescribed by ROTO. 

Strategic head up awareness of ROTO exit turn parameters is provided by an Exit Turn 
symbol shown in true-perspective graphical form, super-imposed on the real world exitway 
visible through the HUD. On final approach, and before ROTO is armed, this symbol is 
shown as a dashed line conformal with the selected ROTO exitway centerline. Upon ROTO 
being armed, this line changes from dashed to solid in form. After touch down, the Exit Turn 
Indication "line" must be modified for clear, unambiguous visibility because of the drastically 
minimized grazing angle resulting from the flight deck's relatively low height above the 
runway. This is accomplished by giving the line a dimension of height (forming a low fence- 
like symbol). Showing the shape of the turn is accomplished by making the 'fence' uniformly 
segmented; at those points where the curve is increasingly more severe in angle (and is also 
relatively farther away from the pilot's eye point) the segments (and gaps between them) will 
be perceived as progressively wider and diminishing in height. (It should be noted here that 
this solution is viewed as provisional, and is not considered completely satisfactory in its 
present version. Chief among the concerns regarding this format symbol are the possibility 
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that this element will clutter the HUD's critically important central display field (even with 
an algorithm for blanking out portions of the element when it competes with other 
symbology), and that it clearly violates the design decision to have dashed lines indicate 
Ready status and solid lines indicated Armed.) This centerline information is supplemented 
by Runway Edge markings to assist the pilot in perceiving position, altitude, track accuracy, 
speed, and perspective cues. 

As an additional source of tactical awareness, Tire Position Indication information is provided 
at the bottom of the HUD. This data indicates current gear position relative to the runway or 
exitway edge and the commanded centerline tracking. This information was provided as a 
safety enhancement employed in correcting for significant deviations from intended centerline 
tracking that could occur in low visibility, with unanticipated runway/exitway conditions, or 
as a result of control errors. 

The final significant modification to a standard HUD format suite is the automatic removal of 
all format elements having exclusively airborne operational functions (e.g., the roll indicator). 
This declutter mode will enable the pilot to more easily attend to and track ROTO's precision 
deceleration, steering, position, and timing. This decluttering function would itself be 
automatically reverted to full-up flight information immediately upon the selection of a go- 
around maneuver. 


Crew Procedures for ROTO Operation 


In this section, a detailed procedural account is provided depicting crew actions and ROTO 
system functioning specifically related to the control and display interfaces for ROTO. This 
description is designed to complement the overall functional timeline developed in the 
operational/mission requirements activity discussed previously (see Appendix A), and to 
better articulate the principal crew procedures and coordination efforts involved in ROTO 
operation. 


A Representative Implementation of ROTO Operation 


On approach into Los Angeles International Airport (LAX), ATC clears the aircraft for 
the ILS (or DGPS) approach to Runway 24 Right. The crew is told to expect to be 
cleared for ROTO later in the approach. The Pilot Not Flying (PNF) accesses the FMS 
page associated with runway ILS selection, specifies the cleared approach information, 
and sets the ROTO system for calculation of appropriate deceleration and steering 
profiles predicated on the system-selected high speed exit. Upon selection of Exit 1, 
the ROTO system enunciates airport, runway, and exitway designations on the 
glareshield (for viewing by both crew members) and on the HUD (and on the head down 
and FMS displays as well). Additionally, ROTO's operational status, "Ready," is 
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indicated. Figure 7.5 portrays this status of the ROTO system while the aircraft is 
turning to intercept the Localizer signal. 

After capture of the localizer and glideslope signals, the aircraft is stabilized on its 
approach to Runway 24 Right, tracking the Glideslope signal in to the predicted touch 
down point. As is apparent in figure 7.6, ROTO is still in Ready status, and the 
runway edge markings and exitway centerline are clearly discernible. The crew awaits 
clearance for ROTO. 

Just prior to the Pilot Flying (PF)'s execution of the flare maneuver in the final seconds 
of flight, the crew has been cleared for the ROTO procedure and the system has been 
armed. Figure 7.7 depicts this change in ROTO system status - ROTO has been armed 
for Automatic (i.e., system-controlled) operation upon touch down, and the exitway 
centerline symbology has consequently turned from dashed to solid line format. The 
PF flares the aircraft and lands. 

Upon touch down, the HUD automatically declutters, leaving symbology relevant only 
to ground maneuvering and to monitoring ROTO progress (or tracking ROTO 
guidance). As is shown in figure 7.8, the exitway centerline symbol has transitioned to 
the fence" format to provide the crew with compelling perspective and angular cues for 
accurately perceiving (and monitoring or tracking) the complex exitway ground path 
dictated by ROTO. Just prior to initiation of the exitway turn maneuver, a positive 
indication of this maneuver s imminent execution is enunciated by a computer-generated 
voice message. Landing gear position is displayed, indicating current position with 
regard to the ROTO-specified centerline of track, and the runway and exitway edges. 

After the aircraft has cleared the runway and continues the commanded deceleration and 
steering profiles on the exitway, the PNF requests clearance to enter the upcoming 
taxiway. In those cases when trailing aircraft are also executing ROTO landings and 
exits, the clearance will need to be timely. In effect, the ROTO procedure is not 
completed until the aircraft has entered the taxiway within the time window required by 
ROTO. 
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8.0 SIMULATED ROTO PERFORMANCE RESULTS 


TIME HISTORIES 


Figures 8.1 - 8. 12 document ROTO time histories using 30 degree spiral exits (reference 2) 
while varying aircraft type, surface condition and crosswind. Each aircraft landed at the mean 
of their dispersion and landing weight as follows: 



MD-11 

MD-81 

Weight (klb.) 

410 

105 

% CG 

23 

16.5 

Airspeed (knots) 

148 

126.5 

Tailwind (knots) (negative is headwind) 

-7.5 

-7.5 

Touchdown Location relative to Threshold (feet) 

1375 

1375 


Table 1 lists the runway occupancy times for these runs. These runs show satisfactory 
braking and steering performance. Runway occupancy time is measured from the runway 
threshold until the aircraft wing tip clears the near side of runway. 

Table 1 


Run # 

Aircraft 

Type 

Surface 

Condition 

Crosswind 

(knots) 

Runway Exit 
Location (feet) 

Runway Occupancy 
Time (sec) 

1 

MD-11 

Dry 

0 

4950 

38.1 

2 

MD-11 

Wet 

15 

6750 

49.5 

3 

MD-11 

Wet 

12.5+/- 2.5 

6750 

49.5 

4 

MD-81 

Dry 

0 

4950 

41.6 

5 

MD-81 | 

Wet 

15 

4950 

41.6 

6 

MD-81 

Wet 

12.5 +/- 2.5 

4950 

41.1 


Note: Positive crosswind direction is from left to right for landing aircraft. A crosswind of increasing 
magnitude causes greater lateral centerline deviation. The simulation studies found a positive 
crosswind caused greater deviation than a negative crosswind. 


Note: Runway Exit Location is relative to runway threshold. 

Note: Runway Occupancy Time is calculated from the time the aircraft crosses the runway threshold until 
the aircraft wing tip clears the near side of the runway. 
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Each simulation run is documented with two pages of time histories. When a plot shares 
more than one variable, the second variable is usually plotted on the right hand Y axis. The 
zero origin of the left and right axis are usually offset so that the variable time histories do not 
cross each other. The X axis of all plots is the runway longitudinal axis in feet. 0 feet is at 
the runway threshold. 


Page 1 ; Bottom Plot 

This plot shows two views of the aircraft position relative to the runway with a right hand 
ROTO turnoff. The left axis shows the aircraft Y position in feet. The runway centerline is 
along the top of the plot. The desired path (dashed line) is along the centerline and then 
curves to the right as the right-handed ROTO exit. Any small perturbations in the dashed 
curves represent exit entrances which the aircraft did not enter. The solid line represents the 
aircraft position. For MD-81 and MD-1 1 simulation runs the first ROTO exit is at position 
3300 feet and 4950 feet, respectively. 

The right axis shows the aircraft Y lateral displacement (solid line) in feet from the runway 
centerline and exit path. The straight-lined funnel shape represents the allowable lateral 
width in which the aircraft can move without running off the pavement. The funnel width is 
the runway and ROTO exit widths minus the aircraft main gear offset, which varies with 
aircraft type. 

Page 1; 2nd from Bottom Plot 

The left axis plots the aircraft ground speed in knots (decreasing trace). The right axis plots 
the aircraft runway occupancy time in seconds. The runway occupancy time at touchdown 
begins at a value greater than zero because it begins counting at the runway threshold. The 

runway occupancy time stops increasing when the aircraft wing tip clears the near side of the 
runway. 


Page 1 ; Middle Plot 

The left axis plots the aircraft lateral acceleration in G’s (lower trace). The right axis plots the 
aircraft lateral jerk in G/sec. Gust cases do not plot the lateral jerk because it is too excessive. 
This study did not ascertain the cause of the gust related jerk (simulation model, control laws, 
sensors) or find a solution for this occurrence. 


Page 1; 2nd from Top Plot 

The left axis plots the aircraft longitudinal acceleration in G’s (lower trace). The right axis 
plots the aircraft longitudinal jerk in G/sec. Gust cases do not plot the longitudinal jerk 
because it is too excessive. This study did not ascertain the cause of the gust related jerk 
(simulation model, control laws, sensors) or find a solution for this occurrence. 
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Page 1; Top Plot 

The left axis plots the percent of main gear brake supply pressure commanded (lower trace). 
When the plot shows 100%, the deceleration command is commanding all of the brake supply 
pressure. The percent of brake supply pressure commanded does not reflect the amount of 
brake supply pressure in use if anti-skid (required by ROTO) is active. 

Please refer to the 3rd and 4th plots on the bottom of plot page 2 for the amounts of available 
mu being used by the main gear. When runway surface friction decreases below that required 
(resulting in skidding), anti-skid decreases brake pressure used just until skidding is alleviated. 
One would not expect 100% supply pressure in use when braking at high speeds on a wet 
surface. The ROTO simulation used in this study implemented the anti-skid function in the 
drag code (for modeling complexity reasons), after its proper location in the brake pressure 
code. 

The right axis plots the aircraft total thrust in pounds (upper trace). 


Page 2; Bottom Plot 

The left axis plots the aircraft rudder position in degrees (lower trace). The right axis plots 
the nose gear position in degrees. 

Page 2; 2nd from Bottom Plot 

The left axis plots the amount of p being used by the aircraft nose gear (lower trace). The 
right axis plots the available aircraft nose gear p.. 


Page 2; 3rd from Bottom Plot 

The left axis plots the amount of p being used by the aircraft main right gear p (lower trace). 
The right axis plots the available aircraft main right gear p. 

Page 2; 4th from Bottom Plot 

The left axis plots the amount of p being used by the aircraft main center gear p (lower trace). 
The right axis plots the available aircraft main center gear |i. 

Page 2; 3rd from Top Plot 

The left axis plots the aircraft track angle relative to the aircraft heading in degrees. The right 
axis plots the aircraft elevator angle in degrees (gradually rising trace). 
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Page 2; 2nd from Top Plot 

The left axis plots the steady tailwind in knots. A headwind would have a negative value. 
The right axis plots the crosswind in knots. If the crosswind is steady it will have a straight 
line value. Gust cases will show a varying crosswind. A positive crosswind blows in a 
negative Y to positive Y direction (left to right as viewed by a landing aircraft). 


Page 2; Top Plot 

The left axis plots the navigation X position data noise content (lower trace). The right axis 
plots the navigation Y position data noise content. 


DISPERSION RESULTS 


The ROTO FORTRAN simulation was modified with an outer loop to allow multiple runs 
Approximately 900 AUTO ROTO runs were simulated that covered the range of expected 
aircraft ground speeds and touchdown locations, spaced 2 knots and 100 feet apart 
respectively. The two data values used for dispersion analysis from these runs were the 
runway occupancy time (ROT) and the ROTO exit location used by the aircraft. 

An effort was then made to calculate the relative probability of one run occurring relative to 
the others. Assuming that the aircraft landing ground speed and touchdown locations are 
normally distributed and independent of each other, the calculations were done as follows: 

1 The aircraft landing ground speed mean and standard deviation were created by adding 
an aircraft’s landing airspeed and expected wind means and variances, respectively, and 
then taking the square root of the summed variance to obtain the standard deviation. 

2. The combined effect of aircraft landing ground speed and touchdown location on the 
relative probability of an individual simulation run occurring was calculated by creating 
a normal probability density function (PDF) for each of the two individual variables. 
This was done by subtracting a cumulative density function (CDF) from the next CDF 
spaced 2 knots apart for the ground speed variable and 100 feet apart for the 
touchdown location variable. A joint PDF, based on the two variables, was created by 
multiplying the individual PDF values together at the intersection values of ground 
speed and touchdown location for each run. A normal CDF is calculated as follows: 

(*-/*) 

1 2<t 2 , 

- 7 — exp v J 

~j27i<y 

The joint normal PDF then represented the relative probability of a run occurring based on 
the aircraft landing ground speed and touchdown location. 

By multiplying a simulation run’s ROT and ROTO exit location by the probability of the run 
occurring, ROT and ROTO exit location statistics could be calculated. Mean, standard 
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deviation and PDF’s were then calculated for ROT and ROTO exit location. A spreadsheet 
was used to manipulate the data and create the figures. 


Figure Descriptions 


There are five types of figures referred to below. The 3-D plots show the ROT value for the 
simulations (e.g. fig. 8-14), the exit location used (e.g. fig. 8-15), and the joint probability 
density function (PDF) of touchdown location and TD ground speed (e.g. fig. 8-16). The 
ROT figure generally shows that the slow and early, and last and fast, landing aircraft have 
the high ROT values. The Exit location figure shows the exit locations used for the 
simulations runs. 

The other two figures are the PDF’s for Exit location (e.g. fig. 8-17) and ROT (e.g. fig. 8-18). 
The Y axis represents the relative occurrences of one outcome (exit used or ROT value) to the 
others. The X axis of the Exit location figure is the possible runway exit locations in feet. 

The X axis of the ROT figure is the possible rounded (to nearest second) ROT values in 
seconds. PDF figures were also created for Exit location and ROT assuming equal 
probabilities for each simulation run. 


Results 


The following table lists the figure numbers for the figures described above: 


Figure Descriptions 

MD-11 

Wet 

MD-11 

Dry 

MD-81 

Wet 

ROT 

8-13 

8-21 

8-27 

Exit Location 

8-15 

8-22 

8-28 

Joint PDF of simulation inputs 

8-16 

8-16 

8-29 

Exit Location PDF 

8-17 

8-23 

8-30 

ROT PDF 

8-18 

8-24 

8-31 

Exit Location PDF assuming equal probabilities 

8-19 

8-25 

8-32 

ROT PDF assuming equal probabilities 

8-20 

8-26 

8-33 

ROT showing incorrect exit predictions 

8-14 




Figure 8-13 uses the same data as figure 8-14, but does not include the spiked ROTs caused 
by 4 incorrect exit predictions. The exit prediction logic used for all the simulation runs was 
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continuously improved during this study, but needs further development to further minimize 
incorrect exit predictions. 

Comparison of the Exit & ROT results show that mean and standard deviations decrease in 
the order of MD-1 1/Wet, MD-1 1/Dry to MD-81/Wet as shown in the table below. The 
lower standard deviation of MD-81 ROT and Exit # taken may be explained by the fact that 
the MD-81 has smaller standard deviations for its touchdown location and TD ground speed. 


Aircraft/Conditions 


Exit # Taken; Exit 0 = 3300 ft ROT (seconds! 
Mean STDDEV Mean STDDEV 


MD-1 1 /Wet 
MD-1 1/Dry 
MD-81 /Wet 


1.65 

0.53 

45.0 

4.04 

1.30 

0.46 

41.3 

3.99 

0.95 

0.40 

41.3 

3.50 


The ROT PDF shows grouped ROTs for each exit location, represented by each line style. 
The left most line style is the 3300 and 4950 foot exit locations, respectively, for the MD-1 1 
and MD-81 aircraft. The Exit location and ROT PDF’s both show that the MD-1 1/Dry 
aircraft make more use of the 4950 foot exit than the MD-1 1/Wet. The MD-8 1/Dry aircraft 
would most likely use the 3300 foot exit more often than the MD-8 1/Wet. The Exit location 

and ROT PDF s using equal probability for all runs show all of the exits being more equally 
utilized. 

The following table shows the percent of MD-1 1 & MD-81 AUTO ROTO aircraft using an 
exit. 


Aircraft/Conditions 

3300 

MD-1 1 /Wet 

0% 

MD-1 1 /Dry 

0% 

MD-81 /Wet 

10.4% 


Exit Locations Used by Aircraft ( feet) 


4950 

6750 

8000 

37.8% 

59.8% 

2.4% 

70.0% 

30.0% 

0% 

84.0% 

5.6% 

0% 


It appears that a 3 ROTO exit (@ 3300, 4950 & 6750 feet past runway threshold) runway 
would satisfy 99.8% (3 sigma dispersion) of the MD-1 1 Diy and MD-81 Wet landings. 
These exits would satisfy ~97.5% (> 2 sigma dispersion = 95.5%) of the MD-1 1 Wet 
landings. 
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9.0 CONCLUSIONS AND RECOMMENDATIONS 


A ROTO architecture and control laws have been designed for a research auto and manual 
ROTO system capable of CATIIIB conditions. System additions to most aircraft would 
include DGPS, FMS exit geometry data, data link, one HUD display, glareshield button(s), 
backdriven tiller, ROTO auto & manual control laws and auto variable braking. To achieve 
the runway occupancy time performance documented in this report; aircraft would also need 
auto-reverse thrust (includes backdriven throttle levers), exit prediction capability and 
consistent reverse thrust stowing techniques based on ground speed. 

Control law designs developed in this study were validated with a 3-DOF aircraft dynamic 
simulation using representative landing gear, engine, rudder actuator, and nose wheel actuator 
models. Approximately 900 simulation runs were performed varying the touchdown point, 
landing speed, aircraft weight and center-of-gravity, wind speeds (headwind/tailwind, 
crosswinds, and gusts), and runway surface conditions (wet & dry). In addition, the 
simulations were run using two different weight class aircraft - an MD-8 1 and MD- 1 1 
aircraft - using the same guidance and control system. This report documents representative 
time histories of these simulations. The study developed ROTO statistics from the 
simulated data including probabilities on runway occupancy time and the use of the given 
available high-speed exits. 

This study recommended ROTO exit positions (@ 3300, 4950 & 6750 feet past runway 
threshold) for MD- 1 1 and MD-8 1 type aircraft used in this study. These exits satisfy from 
2 to 3 sigma touchdown dispersions for aircraft landing weights ranging from light MD-81’s 
up to heavy MD-1 1 ’s. The maximum runway occupancy time for the 6750 foot exit was 
under 53 seconds for a heavy MD-1 1 on a wet/non-grooved concrete surface. This study did 
not investigate other operational airport factors which may not allow for consistent runway 
occupancy times less than 53 seconds. 

This study assumed that ROTO exit entrances would be spiral in nature as described in 
reference 2. This allows for a more smooth transition onto the exit, but adds 1 to 2 seconds 
to the runway occupancy time. It is recommended that ROTO exits be grooved with a 
diamond pattern. If possible, ROTO runways should also be grooved. 

Auto asymmetric braking and fly-by-wire nose gear are no longer required for aft CG MD- 1 1 
aircraft on a wet surface as mentioned in reference 1 . An anti-hysteresis tiller control law, 
commanding the tiller to augment the auto ROTO nose gear, was needed to decrease the 
hysteresis to 2 degrees. It is believed that the tiller control law would cost less than the fly- 
by-wire nose gear alternative. 

Future study is required to find the sensitivity of sub-optimal ROTO designs and operations 
on runway occupancy time. Dynamic piloted evaluation is required of display 
format/guidance elements, look-ahead limits for CATIIIB operations and continuos ROTO 
operations with runway occupancy times less than 53 seconds. Functions required by the 
exit prediction logic need further development: accurate runway friction measurement, 
predicted runway touchdown location/ground speed and general exit prediction algorithm. 
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BASELINE AUTO HIGH-SPEED ROLLOUT / TURNOFF CONTROL SYSTEM 
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The ROTO Computer may he part of the FCC Computer 

The Pilot may ovenide the brake pedals, throttle levers, rudder pedals and/or tiller. 





BASELINE FULL MA NUAL HIGH-SPEED ROLLOUT / TURNOFF CONTROL SYSTEM 



Figure 3.2 
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The ROTO Computer may be part of FCC Computer. 













FIGURE 3.3. ROTO SYSTEM RELIABILITY BLOCK DIAGRAM 
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Figure 16.- Variation ol 990 aircraft braking friction coefficient with ground speed on wet 
grooved and ungrooved surfaces. 5-groove, type VIII, 41 x 15.0-18 main tires: inflation 
pressure, 160 Ib/W. 
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Figure 17.- Variation of 990 aircraft braking friction coefficient with ground speed on 
flooded grooved and ungrooved surfaces. 5-groove, type VIII, 41 x 15.0-18 main tires; 
inflation pressure. 160 Ib/W. 
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Figure 6.2 
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Figure 18.- Variation of 990 aircraft braking friction coefficient with ground speed on wet 
grooved and ungrooved surfaces. Smooth , type VUI, 41 x 15.0-18 main tires; inflation 
pressure, 160 Ib/W. ““ 
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Figure 19.- Variation of 990 aircraft braking friction coefficient with ground speed on 
flooded grooved and ungrooved surfaces. Smooth , type VIII. 41 x 15.0-18 main tires; 
inflation pressure, 160 Ib/in*. _ 
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Figure 6.3 
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ROTO MODE DIAGRAM 
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Figure 7.3. ROTO control and display elements 
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Figure 7.4. ROTO-operational Head Up Display elements 
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Figure 7.5. Terminal area operational elements for 
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OTO implementation: Localizer intercept 




Figure7.6. Terminal area operational elements for ROTO implementation: ILS capture, on glideslope 
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Figure 7.7. Terminal area operational elements for ROTO implementation: At aircraft flare maneuver 
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Figure7.8. Terminal area operational elements for ROTO implementation: Runway roll out 
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Open Integrator 

If (SSE > Max Allowed Deceleration) (9 ft/s/s) 
Reset Integrator to 0 if (SSE < 0) 
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Figure 7.9 Simplified ROTO Auto-Braking Control Law 
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Figure 7.11 Hold Lines 


A/C Y POS. (FT) GROUND SPD (KT) ^ ^ LONG. ACCEL (G) PERCENT BRAKE 



LAT. JRK (GDOT) LONG. JRK (GDOT) 




RUN 1. MD-11 CATIIIB AUTO ROTO (PG 2 OF 2) 

30 DEG SPIRAL EXIT.MID DISP.410KLB.23%CC.-7.5 KNOT TAILWIND 



Figure 8.2 
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A/CYPOS. (FT) GROUND SPD (KT) LAT - ACCEL LONG. ACCEL (G) percent BRAKE 



RUN 2. MD-11 CATIIIB AUTO ROTO (PG 1 OF 2) 

30 DEG SPIRAL EXIT, MID DISP,410KLB.23%CG.WET.15 KNOT CRS WIND 
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Figure 8.3 
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RUN 2. MD-11 CATIIIB AUTO ROTO (PG 2 OF 2) 

30 DEG SPIRAL EXIT.MID DISP,410KLB,23%CG,-7.5 KNOT TAILWIND 
AUTOREVERSE THRUST, ST0WO70 KTS GNO, 4950/6750/8000 EXITS 
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Figure 8.4 
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A/C Y POS. (FT) GROUND SPD (KT) ^ ^ LONG. ACCEL (G) PERCENT BRAKE 



RUN 3. MD-11 CATIIIB AUTO ROTO (PG 1 OF 2) 

30 DEG SPIRAL EXIT.MID DISP.410KLB.23J5CG.WET.12.5+/-2.5 KT CRS WIND 
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Figure 8.5 
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RUN 3. MD-11 CATIIIB AUTO ROTO (PG 2 OF 2) 

30 DEG SPIRAL EXIT, MID DISP.410KLB.23%CG.-7.5 KNOT TAILWIND 
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Figure 8.6 
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A/C Y POS. (FT) GROUND SPD (KT) ^ ^ LONG. ACCEL (G) PERCENT BRAKE 



LAT. JRK (GDOT) LONG. JRK (GDOT) 
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RUN 4. MD-81 CATIIIB AUTO ROTO (PG 2 OF 2) 

30 DEG SPIRAL EXIT.MID DISP.105KLB.175CCG.-7.5 KNOT TAILWIND 



Figure 8.8 
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A/C Y POS. (FT) GROUND SPD (KT) ^ ^ LONG. ACCEL (G) percent BRAKE 
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LAT. JRK (GDOT) LONG. JRK (GDOT) 
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RUN 5. MD-81 CATl 

30 DEG SPIRAL EXIT, MID DISF 
AUTOREVERSE THRUST, STO* 



3 

Ld 

LD 

O 


0.75 

0.50 

0.25 

0.00 

-0.25 

-0.50 

-0.75 



B 

mm 


■HWMS 

BBS 


■ 

mm 

■ B 

MB 

BBI 

■■ 

mim 

fysg 

B 


Bpi 


Bl 



B 

B 

■■ 

■ 

b 

B 



■ 


mm 

Mi 

■ 



O 69- 

Ld 

O 45 J 



[== 




r 












Ld ^ 

o 0 . 







Q 0 

^ -23- 

r 

\ i r \ , 

r\r\ 

= 

n a ^ a 





0 1000 2000 3000 4000 5000 61 


A/C X P0 


Figure 8 


8 



10 


-MAINl MU MAX C MU 




A/C Y POS. (FT) GROUND SPD (KT) ^ ^ LONG. ACCEL (G) percent gpy^ 



RUN 6. MD-81 CATIIIB AUTO ROTO (PG 1 OF 2) 

30 OEG SPIRAL EXIT.MID DISP,105KlB,1755CG.WET,12.5+/-2.5 KT CRS WIND 
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Weight=340K+(480K-340K)*(VEAS-130)/36 
CG=0. 1 2+(0. 34-0.1 2)*(VEAS-130)/36 
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Figure 8.13 
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TD Ground Speed = Airspeed + 
Average Tailwind 


Predict exit prior to TD 

Weight=340K+(480K-340K)*(VEAS-130)/36 
CG=0.12+(0.34-0.12)*(VEAS-1 30)/36 
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Figure 8.14 
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TD Ground Speed = Airspeed 
Average Tailwind 















MD-1 1 Joint Probability Density Function for 2 T.D. inputs (assume independent, gaussian distributed) 
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Figure 8.16 
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MD-11 ROTO Exit Probability Density Function 
Mean=1.65, STDEV=0.53 

W£~r 



MD-11 ROTO Exit Taken 



MD-11 ROTO ROT Probability Density Function 
Mean=45, STDEV=4.04 
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Figure 8.18 
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MD-1 1 Runway Occupancy Time (ROT) seconds 
Curves Represent Exits at 4950, 6750 & 8000 feet 



MD-11 ROTO Exit Probability Density Function 
All landings have equal probability 
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Figure 8.19 
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MD-11 ROTO Exit Taken 



MD-1 1 ROTO ROT Probability Density Function 
Ail landings have equal probability 



MD-11 Runway Occupancy Time (ROT) seconds 
Curves Represent Exits at 4950, 6750, 8000 & 10000 feet 
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Figure 8.21 
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MD-11 ROTO Exit Probability Density Function 
Mean=1.3, STDEV=0.46 
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Figure 8.23 
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MD-11 ROTO Exit Taken 
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MD-11 Runway Occupancy Time (ROT) seconds 
Curves Represent Exits at 4950, 6750 & 8000 feet 
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MD-11 ROTO Exit Taken 



MD-1 1 ROTO ROT Probability Density Function 
1 20 | All landings have equal probability 
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Figure 8.26 
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MD-1 1 Runway Occupancy Time (ROT) seconds 
Curves Represent Exits at 4950, 6750, 8000 & 10000 feet 
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Predict exit prior to TD MD-81 ROTO Exit Used Wet, Exits=3300,4950,6750, 8000 

Autoreverse Thrust 

Weight=82K+(128K-82K)*(VEAS-1 10)/33 ! Stow Reverse Thrust=70 kt gd 

CG=-0.008+(0.334-(-0.008))*(VEAS-1 10)/33 ^ 












MD-81 Joint Probability Density Function for 2 T.D. inputs (assume independent, gaussian distributed) 
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MD-81 ROTO Exit Probability Density Function 
Mean-0.95, STDEV-0.4 
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Figure 8.30 
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MD-81 ROTO Exit Taken 



MD-81 ROTO ROT Probability Density Function 
Mean=41.3, STDEV=3.50 
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Figure 8.31 


MD-81 Runway Occupancy Time (ROT) seconds 
Curves Represent Exits at 3300, 4950 & 6750 feet 



MD-81 ROTO Exit Probability Density Function 
All landings have equal probability 
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MD-81 ROTO Exit Taken 



MD-81 ROTO ROT Probability Density Function 
All landings have equal probability 



MD-81 Runway Occupancy Time (ROT) seconds 
Curves Represent Exits at 3300, 4950, 6750 & 8000 feet 




Table 3.1. SUMMARY - HAZARD CLASS I 


SYSTEM FUNCTION 
Runway Acquisition Clearance 

ROTO Steering Guidance 
ROTO Steering Guidance 
ROTO Steering Guidance 


HAZARD DESCRIPTION 
Object on runway 

Erroneous steering guidance 
Erroneous rudder 
Erroneous nosegear steering 


FAILURE CONDITION 
Possible crash 

Loss of control 
Loss of control 
Loss of control 



Table 3.2. SUMMARY - HAZARD CLASS II 


SYSTEM FUNCTION 

HAZARD DESCRIPTION 

FAILURE CONDITION 

ROTO Navigation 

Loss of navigation data 
(GPS, IRS, airport data base) 

Loss of ROTO capability 

Loss of control if loss 
occurs at critical time 

ROTO Navigation 

Erroneous navigation data 

Loss of control 

ROTO Steering Guidance 

Loss of all guidance command 

Loss of ROTO capability 

Loss of control if loss 
occurs at critical time 

ROTO Steering Guidance 

Loss of rudder 

Loss of ROTO capability 

Loss of control if loss 
occurs at critical time 

ROTO Steering Guidance 

Loss of nosegear 

Loss of ROTO capability 

Loss of control if loss 
occurs at critical time 

ROTO Deceleration 

Loss of deceleration command 

Loss of ROTO capability 

Loss of control if loss 
occurs at critical time 

ROTO Deceleration 

Erroneous deceleration 
commands 

Loss of ROTO capability 

Loss of control if loss 
occurs at critical time 

ROTO Deceleration 

Erroneous auto brake 

Loss of ROTO capability 

Loss of control if loss 
occurs at critical time 

ROTO Deceleration 

Erroneous reverse thrust 

Loss of ROTO capability 

Loss of control if loss 
occurs at critical time 

ROTO Display 

Loss of pilot display 

Loss of ROTO capability 
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occurs at critical time 
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Loss of aircraft weight 

Loss of ROTO capability 
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Table 3 . 3 . SUMMARY - SYSTEM ARCHITECTURE RECOMMENDATIONS/CONSIDERATIONS 
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Table 3.4. SUMMARY - SYSTEM ARCHITECTURE RECOMMENDATIONS/CONSIDERATIONS 

HAZARD CLASS 1 1 
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Failure The ROTO system failure monitoring system must be designed so that no single 
monitoring failure or combination of failures not shown to be extremely improbable will cause 

system loss of the failure monitoring system (undetected failure, misleading information, 



Table 3.4. SUMMARY - SYSTEM ARCHITECTURE RECOMMENDATIONS/CONSIDERATIONS 

HAZARD CLASS 1 1 (continued) 
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MAKEFILE is a file used by MAKE to create C object files required by 
MODELROTOBIG.M, which are then converted to CMEX files. 

CLIBS = -lm -lmalloc -lc 

CFLAGS = +z -z -Aa -Ae -g -DDAYFMEX -I$(HOME)/matlab/extem/include 
atmos : funcnosetireforces.o 

cmex funcnosetireforces.o 
.f.o: 

cc S(CFLAGS) funcnosetireforces.c 

atmos : funcinput.o 

cmex funcinput.o 
.f.o: 

cc $(CFLAGS) funcinput.c 

atmos : funcinputnav.o 

cmex funcinputnav.o 

.f.o: 

cc $(CFLAGS) fiincinputnav.c 
atmos : funcinterpr.o 

cmex funcinterpr.o 
.f.o: 

cc $(CFLAGS) funcinterpr.c 
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RUNROTO.M is the main MATLAB script file to run the ROTO simulation from the 
MATLAB command line. It calls the MATLAB SIMULINK model 
MODELROTOBIG.M, which you may view by typing modelrotobig from the 
MATLAB command line. The following pages contain the input MATLAB script files 
and exit prediction algorithm files. 

% Initialize all input variables 
clear all; 
format short e; 

% if you were to do multiple runs, do them here with clear variables; 

% this will not reload the model and functions 

centerline; 

centerline22; 

% generic airplane 
inputdef; 

% unique mdl 1 characteristics 
indefmdll; 

% run specifics 
inmdl 170; 

% some additional initializations 
beforeloop; 
save mdl 170; 
createarrays; 
globalfile; 

% Predict exit selection before touchdown 

% create U_TD,UT _TD,X_TD; input files currently set these values, 

% In reality they would need to be predicted on the glideslope 
UJTD=U; 

UT_TD=UT; 

X_TD=X; 

VWSS_TD=VWSS; 

lastnexit=0; 

NONROTO-O; 

while ((lastnexit~NEXIT)&(NEXIT<5)) 
lastnexit=NEXIT; 

aborteariy=funcpredictabort(U_TD,UT_TD,X_TD,VWSS_TD,VEXIT,EXITPOS,NEXIT); 

if(abortearly) 

NEXIT=NEXIT+1; 

end; 

if(NEXIT>4) 

NONROTO=1; 

end; 

EXITPOS=EXITLOC(NEXIT); 

end; 


% ROTO Simulation 
ifi[NONROTO=0) 

% Run roto simulation 

T1ME=0; 

timt=0; 

TSTEP=0; 

[T,AEROX,AEROY]=Iinsim('modelrotobig',70,[],tDT,DT,DT,0,3,2]); 

% Write out results 
fid=fopen('nose.output 1 .data’, V); 
scriptwriteplotdata; 
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else 

dispCNON-ROTO LANDING'); end; 



VARIABLE DEFINITIONS 


Variable definitions used in the ROTO simulation. 


The following defines output variables: 


Variable 

X 

Y 

YSTRIPE 

U 

PS1 

DV 

DU 

DY 

DPSI 

SIGSTR 

SIG 

DELRC 

DPERP 

VPERP 

DDPSI 

FNOSE 

FC 

FM 

PMM 

PSIP 

THST, THRUST 
TR 

MNAERO 

PSIC 

NC 

PSICGV 

MUMAX 

MUMAXR 

MUMAXL 

MUMAXC 

MUNOSE 

MUR 

MUL 

MUC 

LATG 

LONGG 

LATJ 

LONGJ 

YAW 

DRU 

DRUR 

DRUDEAD 

DRMPUU 

DECELCOM 

NEXIT 

RUNWAYTIME 
POSEXIT,NEGEXIT 
PERCENT BRK 
TAILWIND 
CROSSWIND 


Definitions 

Aircraft longitudinal position on the runway 
Aircraft lateral position on the runway 
Runway and exit centerline lateral position 
Aircraft ground speed in knots (inside sim it is fps) 

Aircraft heading 
Aircraft lateral accelleration 
Aircraft longitudinal accelleration 
Aircraft lateral deviation rate 
Aircraft heading rate 

Aircraft nose wheel angle relative to aircraft heading 

Aircraft nose wheel angle 

Rudder command 

Aircraft lateral deviation 

Aircraft lateral deviation rate 

Aircraft heading acceleration 

Aircraft nose tire side force 

Aircraft center tire side force, after limited to avail, friction 

Aircraft R&L tire side force, after limitied to avail, friction 

Aircraft R&L tire vertical force 

Aircraft heading relative to centerline heading 

A ircraft total thrust 

Aircraft rudder moment 

Aircraft side slip moment 

Aircraft center gear heading 

Aircraft center gear cornering power 

Aircraft center gear heading relative to aircraft velocity vector 

Runway surface available friction coef. at nose gear 

Runway surface available friction coef. at right main gear 

Runway surface available friction coef. at left main gear 

Runway surface available friction coef. at center main gear 

Friction used by nose gear 

Friction used by right main gear 

Friction used by left main gear 

Friction used by center main gear 

Aircraft lateral G 

Aircraft longitudinal G 

Aircraft lateral jerk 

Aircraft longitudinal jerk 

Aircraft yaw angle 

Aircraft upper rudder angle 

Aircraft upper rudder angle rate 

Aircraft upper rudder value after deadband 

Aircraft upper rudder mod piston position 

Aircraft symmetrical deceleration command 

Designnated exit # (1,2, or 3) 

Aircraft runway occupancy time 

Aircraft main gear must stay within these lateral bounds 

Avg. Percent of maximum brake pressure used by all main gear 

knots 

knots 
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XNAV 

YNAV 

ELEV 


DGPS noise added to actual aircraft X position 
DGPS noise added to actual aircraft Y position 
Aircraft elevator angle 
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The following describes input contants: 
* 


Variable 

* A 

* ASEFF 

* B 

* BC 

* BW 

* CMR 

* C 

* CDRAG 

* CLIFT 

* CMOM 

* CNB 

* CYB 

* DELB 

* DELS 
*DECELSET 

* HCG 

* HS 

* IYAW 

* LMAC 

* LTA1L 

* NM 

* OD 

* RB 

* RP 

* RS 

* SP NLG 

* SPM 

* S 

* SW 

* THETA 

* TREAD 

* VEAS 

* VW 

* VTW 

* W 

* WET 

* ICE 

* SNOW 

* SLUSH 

* FLOOD 

* WS 

* FLATMU 

* KSI 

* KS2 

* KS3 

* PSUP 

* PRET 

* TMGD 

* TNGD 

* TSPOIL 

* TDELAY 
♦RRPHASE1 
♦RRPHASE2 
*KBPHASE1 
♦KBPHASE2 


Definition 

DISTANCE-NLG TO CG 
ANTISKID EFFICIENCY 
DISTANCE-CG TO MLG 
DISTANCE-CG TO CLG 
DISTANCE - CG TO 
RUDDER YAW MOMENT COEFF. 

LIFT MOMENT ARM TO CG 
A/C AERODYNAMIC DRAG COEFF. 

A/C AERODYNAMIC LIFT COEFF. 

A/C AERODYNAMIC MOMENT COEFF. 
AERODYNAMIC SIDE-SLIP MOMENT COEFF. 
AERODYNAMIC SIDE-SLIP FORCE COEFF. 
NLG TIRE DEFL., AT LOAD RB 
NLG TIRE DEFL., RATED 
AUTOBRAKE DECELERATION SETTING 
CG HEIGHT 

NLG TIRE SECTION HEIGHT 

A/C YAW MOMENT OF INERTIA (ABOUT CG) 

LENGTH OF THE MAC 

DISTANCE-CG TO TAIL CP 

MLG TIRE CORNERING POWER (PER TIRE) 

NLG TIRE OD 


FT 

FT 

FT 

FT 

l./DEG 

FT 


IN 

IN 

FT/SEC 

FT 

IN 

SLUG-FT**2 

FT 

FT 

LB/DEG 

IN 


NLG TIRE VERT LOAD AT DEFL. DELB (FOR SPRING RATE CALC.) LB 


NLG TIRE RATED PRESS. (LOADED) PSI 

NLG TIRE RATED LOAD LB 

TIRE STATIC PRESSURE (LOADED) PSI 

MLG TIRE STATIC PRESSURE (LOADED) PSI 

NOSEWHEEL SPACING IN 

A/C WING AREA FT**2 

NLG STRUT CANT ANGLE DEG 

DISTANCE BETWEEN MLG'S FT 

INITIAL AIRSPD A/C VELOCITY KTS 

CROSSWIND VELOCITY KTS 

TAILWIND VELOCITY KTS 

A/C WEIGHT LB 

SWITCH TO USE WET RUNWAY FRICTION (1 IS ON) — 

SWITCH TO USE ICY RUNWAY FRICTION (2 IS ON) 

SWITCH TO USE SNOWY RUNWAY FRICTION (3 IS ON) — 

SWITCH TO USE SLUSHY RUNWAY FRICTION (4 IS ON) 

SWITCH TO USE FLOODED RUNWAY FRICTION (5 IS ON) — 

NLG TIRE SECTION WIDTH IN 

SWITCH FOR FLAT CORNERING MU VS. SKID VEL. CURVE (1 IS FLAT) 

(STEERING VALVE SPOOL DISPLACEMENT)/(STEER ERR) IN/DEG 
(STEERING RATE)/( VALVE FLOW) (DEG.SEC)/(IN**3/SEC) 

(STEERING ACTUATOR PRESSURE)/(STRT GD MOM) (PSI/IN-LB) 

HYDRAULIC SUPPLY PRESSURE PSI 

HYDRAULIC RETURN PRESSURE PSI 

MAIN GEAR TOUCHDOWN TIME SEC 

NOSE GEAR TOUCHDOWN TIME SEC 

TIME BETWEEN MAIN GEAR TOUCHDOWN & SPL DEPLOYMENT SEC 
TIME BETWEEN SPOILER DEPLOYMENT & START OF BRKE RP SEC 
PHASE 1 BRAKE PRESSURE RAMP RATE PSI/SEC 

PHASE 2 BRAKE PRESSURE RAMP RATE PSI/SEC 

PHASE 1 BRAKE PRESSURE GAIN (PSI/SEC)/(FT/SEC**2) 

PHASE 2 BRAKE PRESSURE GAIN (PSI/SEC)/(FT/SEC**2) 
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* NWWLG 

* NWCLG 

* BDEXP 

* BDP 

* BDK 


NUMBER OF WING GEAR WHEELS 

NUMBER OF CENTER GEAR WHEELS 

CONSTANT USED TO DEFINE BRAKE DRAG VS. PRESSURE CURVE - — 
CONSTANT USED TO DEFINE BRAKE DRAG VS. PRESSURE CURVE PSI 
CONSTANT USED TO DEFINE BRAKE DRAG VS. PRESSURE CURVE LBS/PSI 
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CENTERLINE.M is a file that gives coordinates for the spiral exit centerline every 2 feet 
apart along its length and is called by RUNROTO once. Because of it many pages of data 
it is not included in this report. NASA has the file in electronic form. This file is not 
used by the aircraft navigation. It is used in the simulation by FUNCINPUT to find the 
actual centerline lateral error of the aircraft for performance evaluation purposes. 
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CENTERLINE22.M is a file containing the coordinates and radius of the spiral exit 
centerline and is called by RUNROTO once. The points are spaced so that the 
perpendicular error from the actual path is less than 0.5 feet. There are about 26 points. 
The points get closer together as the path radius decreases. The database requirements 
section recommends that the exit’s X,Y coordinates be in units of 1/2 foot using 2 byte 
integers. This particular input file has stored the X,Y coordinates in units of feet, but the 
X,Y coordinates have no greater resolution than 1/2 foot. 

spiralexit22 = [ 

O.OOOOOD+OO... 

O.OOOOOD+OO... 

0.32767D+5... 

0.17600D+03... 

0.5D+00... 

0.10634D+05... 

0.354D+03... 

0.400D+01 ... 

0.52872D+04... 

0.490D+03... 

0.105D+02... 

0.38197D+04... 

’ 0.601 5D+03.. . 

0.195D+02... 

0.31091D+04... 

0.7105D+03... 

0.32D+02... 

0.26287D+04... 

0.8055D+03... 

0.47D+02... 

0.23164D+04... 

0.9075D+03... 

0.675D+02... 

0.20523D+04... 

0.985D+03... 

0.865D+02... 

0.18868D+04... 

0.10695D+04... 

0.1 1 15D+03... 

0.17330D+04... 

0.1 145D+04... 

0.1375D+03... 

0.16135D+04... 
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0.1219D+04... 

0.1675D+03... 

0.15094D+04... 

0.I290D+04... 

0.201D+03... 

0.14201D+04... 

0.1343D+04... 

0.229D+03... 

0.I3582D+04... 

0. 1362D+04... 
0.2395D+03... 
0.I3369D+04... 

0.1362D+04... 

0.2395D+03... 

0.32767D+5... 

0. 1800D+04... 
0.4925D+03... 
0.32767D+5... 

0. 1800D+04... 
0.4925D+03... 
-.80000D+03... 

0.18495D+04... 

0.519D+03... 

-.80000D+03... 

0.1901D+04... 

0.5415D+03... 

-.80000D+03... 

0.19555D+04... 

0.5615D+03... 

-.80000D+03... 

0.2009D+04... 

0.577D+03... 

-.80000D+03... 

0.2064D+04... 

0.588D+03... 

-80000D+03... 

0.21 195D+04... 

0.596D+03... 

-.80000D+03... 

0.21775D+04... 

0.5995D+03... 

-.80000D+03... 



0.21 995 D+04... 
0.60000D+03... 
-.80000D+03... 


0.22015D+04... 

0.60000D+03... 

0.32767D+5... 


* 

]; 

spiralexitx22=spiralexit22(:, I ); 
spiralexity22=spiralexit22(:,2); 
spiralexitr22=spiralexit22(:,3); 
clear spiralexit22; 
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INPUTDEF.M is a MATLAB script file which initializes the default simulation variables 

for all aircraft types and flight conditions and is called by RUNROTO once. 
false=0; 
true=l; 

ASEFF = 0.75; 

DT = 1.0E-03; 

DECELSET = 0.0; 

ENDTIM = 350.0; 

INTVL = 0.5; 

VEAS = 166.; 

VSTOP = 1.0; 

VWSS= 15; 

VTWSS = 10.; 

WET = 1; 

MU ROLL = 0.15; 

PSUP = 3000.; 

PRET = 60; 

TMGD = 0.0; 

TNGD = 6.0; 

TSPOIL = 1.3; 

TDELAY = 3.0; 

RRPHASE1 = 400.; 

RRPHASE2 = 1200.; 

KBPHASE1 = 600.; 

KBPHASE2 = 1800.; 

FLATMU = 1; 

RUNWA YWIDTH= 1 50; 

EX1TWIDTH=90; 

% table lookup uses last four points to extrapolate out of table, if you want 
% a limit in the table, use four points at each end with that limit. 

MUSKIDT = [ 

O.OOOOOOOOOE+OO; 

20 ; 

40; 

60; 

100 ; 

100 . 0 ; 

1 ; 

MUSKID = [ 

1 . 0 ; 

0 . 6 ; 

0.42; 

0.36; 

0.32; 

0.32; 

1 ; 

STEERT = [ 

-23; 

-18.9; 

- 12 . 6 ; 

-6.3; 

-3.15; 

0 . 0 ; 

3.15; 

6.3; 

12 . 6 ; 

18.9; 

23; 


A-12 



]: 

STEER = [ 

8.4; 

7.; 

4.85; 

2.55; 

1.3; 

0 . 0 ; 

-1.3; 

-2.55; 

-4.85; 

-7.; 

-8.4; 

1 ; 

WETMUT = [ 
0 . 0 ; 

20 .; 

40.; 

50.; 

70.; 

80.; 

100 .; 

120 .; 

140.; 

160.; 

180.; 

]; 

WETMU = [ 
0.941; 

0.533; 

0.378; 

0.333; 

0.282; 

0.271; 

0.266; 

0.266; 

0.266; 

0.27; 

0.274; 

]; 

ICEMUT = [ 

0 . 0 ; 

200 .; 

]; 

1CEMU = [ 
0.27E-01; 
0.27E-01; 

]; 

SNOWMUT = [ 

0 . 0 ; 

200 .; 

]; 

SNOWMU = [ 
0.162; 

0.162; 

]; 

SLUSHMUT = [ 
0.0; 
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200 .; 

]; 

SLUSHMU = [ 
0.55; 

0.55; 

]; 

FLOODMUT = [ 

0 . 0 ; 

10 .; 

20 .; 

30.; 

40.; 

50.; 

60.; 

70.; 

80.; 

90.; 

100 .; 

110 .; 

120 .; 

130.; 

140.; 

150.; 

160.; 

170.; 

180.; 

190.; 

200 .; 

]; 

FLOODMU = [ 
0.929; 

0.645; 

0.465; 

0.353; 

0.276; 

0 . 22 ; 

0.18; 

0.147; 

0.119; 

0.98E-01; 

0.84E-01; 

0.72E-01; 

0.62E-0 1 ; 
0.55E-01; 

0.53 E-01; 
0.5E-01; 

0.51 E-01; 
0.52E-01; 
0.53E-01; 
0.54E-01; 
0.55E-01; 

]; 
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INDEFMD1 l.M is a file initializes MD-1 1 aircraft characteristics independent of flight 

contidion and is called by RUNROTO once. 

A = 78.2562256; 

B = 2.45735645; 

BC = 5.00903320; 

BW = 165.369904; 

C = 2.21834278; 

CDRAG = 0.1651; 

CLIFT = 0.29; 

CMOM = 0.160003304E-01; 

CMR = -0.262E-02; 

CNB = 0.37E-02; 

CYB = -0.24E-01; 

DELB - 1.1; 

DELS = 3.4; 

HCG = 2.71658039; 

HS = 9.8; 

IYAW = 26000000.0; 

LMAC = 24.6479797; 

LTAIL = 83.7390137; 

NC = 4426; 

NM = 4806.0; 

OD = 39.6; 

RB = 8000; 

RP = 203; 

RS = 39500; 

SP= 167; 

SPM = 188; 

S = 25.; 

SW = 3647.5; 

W1NGSPAN=83; 

TREAD = 34.677; 

THETA = 9.5; 

W = 480000.; 

KS1 = 0.873E-02; 

KS2 = 0.965; 

KS3 = 0.842E-02; 

BDEXP=0.7; 

BDP=125: 

BDK= 190.6; 

NWWLG = 4; 

NWCLG = 2; 

WS= 15.5; 

% table lookup uses last four points to extrapolate out of table, if you want 
% a limit in the table, use four points at each end with that limit. 

THS1T = [ 

0 . 0 ; 

66.1249695; 

132; 

]; 

THSl = [ 

15048; 

10602; 

7335; 

1; 

THSRT = [ 

85.0; 

123.84; 
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126; 

135; 

140; 

]; 

THSR = [ 

-28740; 

-40700; 

-40900; 

-41460; 

-42630; 

]; 

THSTIT = [ 
9.5; 

7.0; 

5.5; 

3.5; 

0 . 0 ; 

]; 

THSTl = [ 

-39350; 

-23340; 

- 7060; 

6800; 

6750; 

]; 

THSNT = [ 
0 . 0 ; 

66.125; 

132; 

]; 

THSN = [ 
7581; 

4680; 

2535; 

]; 

THSRNT = [ 
0 . 0 ; 

60.889; 

65.65; 

84.96; 

123.84; 

125.97; 

134.85; 

140.18; 

147.29; 

150.0 ; 

]; 

THSRN = [ 
-1400.0; 
-6350.0; 
-6750.0; 
-8425.0; 
-11900.0; 
-12150.0; 
-13050.0; 
-14200.0; 
-15475.0; 
-16200.0; 
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]; 

GAMMAT = [ 
200000.; 
257500.; 
280000.; 
300000.; 
330000.; 
365000.; 
405000.; 
440000.; 
480000.; 
510000.; 
550000.; 

]; 

GAMMA = [ 
0.224; 

0 . 2 ; 

0.192; 

0.186; 

0.178; 

0.17; 

0.162; 

0.156; 

0.15; 

0.146; 

0.142; 

]: 

VALVET = [ 
-0.13; 

- 0 . 12 ; 

- 0 . 11 ; 

- 0 . 11 ; 

-l.E-01; 

-0.9E-0 1 ; 

-0.8E-0 1 ; 

-0.7E-01; 

-0.6E-01; 

-0.5E-01; 

-0.4E-01; 

-0.3E-01; 

-0.2E-01; 

-0.1E-01; 

0 . 0 ; 

0.1E-01; 

0.2E-01; 

0.3E-01; 

0.4E-01; 

0.5E-0I; 

0.6E-01; 

0.7E-01; 

0.8E-01; 

0. 9E-01; 

1. E-01; 

0 . 11 ; 

0 . 11 ; 

0 . 12 ; 

0.13; 

]; 
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VALVE = [ 
-0.91; 

-0.91; 

-0.91; 

-0.91; 

-0.875; 

-0.765; 

- 0 . 6 ; 

-0.455; 
-0.34; 
-0.245; 
-0.165; 
-0.9E-0 1 ; 
-0.3E-01; 
0 . 0 ; 

0 . 0 ; 

0 . 0 ; 

0.3E-01; 

0.9E-01; 

0.165; 

0.245; 

0.34; 

0.455; 

0 . 6 ; 

0.765; 

0.875; 

0.91; 

0.91; 

0.91; 

0.91; 

]; 
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INMD1 170.M is a file initializes variables specific to a flight condition and is called by 

RUNROTO once. 

MDl l=true; 

DT = 0.005; 

DECELSET = 0.0; 

ENDTIM = 100; 

INTVL = 0.5; 

% veas range 130 kt to 166 kt 
VEAS = 166; 

VTWSS = 10.; 

VWSS= 15; 

% weight range 340k to 480k 
W = 480000; 

WET = 1; 

KYAWRT = 7.75; 

KYDI = 5; 

KY= 12; 

KRAD = -40000; 

XSTART = 2500; 

RCMDLM = 1000.0; 

EXSTOP = 1.0; 

RADTAU = 1.0; 

IG3 = 0.1; 

IG5 = 0.002; 

1G16 = 0.06; 

IG14 = 23.0; 

IKBANDRD = 0.85; 

IKMECH = 0; 

NSDEAD = 0.01; 

NSRTLM = 0.3; 

NSPSLM = 8.0; 

KBANDNS = 0.7; 

DECLIM = 9.0; 

COMPUP = 0.05; 

ILATBIAS = 0.0; 


ITURB = false; 

NAVUP = 0.1; 
IRANDM = false; 
ACCURRMS = 4.0; 
ACCURTAU = 30.0; 


% CG range .12 to .34 
CG = 0.34; 

% CALCULATE NEW CG VARIABLES BASED ON FRACTION CG 
CGIN=131 1.947+CG*295.779; 

A=(CGIN-473.437)/12.; 

B=( 1 442-CG IN)/ 1 2. ; 

BC=( 1 472.62-CGIN)/ 1 2.; 

C=(CGIN-(1311.947+0.25*295.779))/12.; 

HCG=(209.32-(sqrt(CGIN A 2+(-2 1 ) A 2))*sin(atan(2 l/CGIN)+0.00 1 93))/12.; 
LTAIL=(2417.38-CGIN)/12.; 

% SET 1H BASED ON CG 
PCTMAC=100.*CG; 

ELEV= 8; 

CGFWD=12; 

CGAFT=34; 

CLFWD=0.123; 

CLAFT=0.226; 

CMFWD=0.515; 
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CMAFT=0.216; 

CDFWD=0.1746; 

CDAFT=0.1651; 

DCLDE=0.008; 

DCMDE=-0.025; 


% other constants 
BUFSIZ = 0.15; 

% DT WOULD HAVE BEEN COMPUP IF IT WERE RUNNING AT COMPUP 
% RADLAG=exp(-DT/RADTAU); 

RUDLIM = 23.; 

C5 = 0.07; 

KBYR2 = 7.75; 

KLAG = 0.9; 

GSPFLL = 130.0; 

R2D = 180./3. 141593; 

REXP2 =exp(-DT/ 10.0); 

REXP3 =exp(-DT/ 1.0); 

Cl 6 = .084; 

RUDRLM = 30.0; 

KBYRDT=10.0; 

VTW=VTWSS; 

AUTOREVERSE=l ; 

% von karman filter 

WNYQ = 3.1416/DT; 

SRNYQ = sqrt(WNYQ); 

WD = 647.9/ 1750.; 

SG=VWSS/1.5; 

GD = SG*sqrt( 1 ,0/( WD*3 . 1416)); 

WF = 100. * WD; 

% VON KARMAN COEFFICIENTS 
VAPB=(1.339+1.1 18)/WD; 

VATB= 1 .339* 1.11 8/WD/WD; 

VCPD=(0. 1 277+0.0 1 46)/WD; 

VCTD=0. 1277*0.0 146/WD/WD; 

VEPF=(2. 1 87+0. 1 833)/WD; 

VETF=2. 1 87*0. 1 833/WD/WD; 

VG=0.021/WD; 
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BEFORELOOP.M is a file initializing additional inputs prior to the simulation beginning 
and is called by RUNROTO once. This file is not aircraft type specific. Some if 

statements are aircraft type specific. 

% FORCE MU TO BE FLAT FOR NON-DRY CONDITIONS 
if(WET>0) 

FLATMU=true; 

else 

FLATMU=false; 

end; 

G3=IG3; 

G20=IG3; 

G5=IG5; 

G16=IG16; 

G14=IG14; 

KBANDRD=IKBANDRD; 

LATBIAS=ILATBIAS; 

VTW=VTWSS; 

VW=VWSS; 

CLIFT=CLFWD+((CLAFT-CLFWD)/(CGAFT-CGFWD))*(PCTMAC-CGFWD)+DCLDE*ELEV; 

CDRAG=CDFWD+((CDAFT-CDFWD)/(CGAFT-CGFWD))*(PCTMAC-CGFWD); 

CMOM =CMFWD+((CMAFT-CMFWD)/(CGAFT-CGFWD))*(PCTMAC- 
CGFWD)+DCMDE*ELEV; 


% INITIALIZATION AND CALCULATIONS 
PI =3.141592654; 

DPSI=0.0; 

PSI= 0.0; 

% FOR SOME REASON, EXITPOS VARIABLE DOES NOT LIKE TO BE NEGATIVE 
% MAKE XSTART SO THAT FIRST EXIT IS AT 0 
X=XSTART; 

Y=0.0; 

SIG=0.0; 

DELRC=0.0; 

Kl= (RS-RB)/(DELS-DELB); 

K2 = RB/DELB; 

SST=S*sin(THETA/57.3); 

D1 = DELB-RB/K1; 

% U IS A GRND SPD 

U=(VEAS+VTWSS)* 1 .689; 

V=0.0; 

DX = U*cos(PSI) - V*sin(PSI); 

DY = V*cos(PSI) + U*sin(PSI); 

TRQLAG=exp(-DT/0.005); 

COUNT=0.0; 

ACCURLAG=exp(-DT/ACCURTAU); 

ASYMBRAKE=false & (NWCLG~=0); 

% XYACCUR SETUP 

ACCURRMS = sqrt(3 .141 6/DT)* ACCURRMS; 

BUFSIZ=0.15*0; 

% patch constants 

WETPATCH= WET ; 

PATCHNOW=false; 

PATCHLEN=10; 

PATCHTYPE=5 ; 

FIRSTPATCH=0; 
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LASTPATCH= 11500; 

PATCHOFTEN=10; 

DISPLN=2250; 
if( MD 1 1 ) 
if(DISPLN<l 150) 

EXITLOC( 1 )=4950; 

EXITLOC(2)=6750; 

EXITLOC(3)=8000; 

EXITLOC(4)= 10000; 
else 

EXITLOC( 1 )=4950; 

EXITLOC(2)=6750; 

EXITLOC(3)=8000; 

EXITLOC(4)= 10000; 
end; 
else 

if(DISPLN<l 150) 

EXITLOC(1)=3300; 

EXITLOC(2)=4950; 

EXITLOC(3)=6750; 

EXITLOC(4)=8000; 

EXITLOC(5)= 10000; 
else 

EXITLOC(1)=3300; 

EXITLOC(2)=4950; 

EXITLOC(3)=6750; 

EX1TLOC(4)=8000; 

EXITLOC(5)= 10000; 
end; 
end; 

%other 

AUREVSEC= 3.0; 

AUREVTHIRD= 1; 
if(ITURB=l) 

RUDDEF= 0; 
else 

RUDDEF=VWSS; 

end; 

DECELON=false; 

LASTDECCR=0; 

LASTDECCL=0; 

DECELSET=0; 

REVERSE=true; 

AUTOREV=true; 

USEMUABORT=true; 

NEXIT=1; 

EXITPOS=EXITLOC(NEXIT); 

% ASSUME PLANE IS DECRABED 

UT =(U - sin(0)*(VWSS* 1.689) - cos(0)*(VTWSS*I.689)V1.689; 
VEXIT= 70.0*1.689; 

INITRUNTIME= 1 /(U/X); 

OUT A RRSIZE=ENDT1M/DT; 

TSTEP=0; 

timt=0; 

TIMEOCCLAG=(9000-X)/((VEXIT+Uy2^>fINITRUNnME; 

S WITCHDI V=0.000000000 1 ; 

ALSAMPLE=DT ; 
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LOWSPD=20* 1.689; 
BRKBUF=100; 
MANUAL=0; 
DRATE=8.0 ; 
DTEMP=7.5; 
DTIME=41; 



CREATEARRAYS.M is a file which declares array variables (mostly output), which will 
decrease the simulation run time. It called by RUNROTO once. 

RWFC = zeros(size( 1:5000)); 
for i = 1:5000 
RWFC(i)=WET; 
end; 

breakoutin=(BDP:3000+BDP); 
breakout=BDK *(breakoutin-BDP). A BDEXP; 

XTH = zeros(OUTARRSIZE, 1 ); 

DXTH = zeros(OUTARRSIZE, 1 ); 

DUTH = zeros(OUTARRS!ZE, 1 ); 

UTH = zeros(OUTARRSIZE, 1 ); 

YTH = zeros(OUTARRSIZEJ); 

DYTH = zeros(OUTARRSIZE, 1); 

DVTH = zeros(OUTARRSIZE, 1 ); 

PSITH = zeros(OUTARRSIZE,l); 

DPSITH = zeros(OUTARRSIZE, 1 ); 

DDPSITH = zeros(OUTARRSIZE, 1 ); 

S1GTH = zeros(OUTARRSIZEJ); 

SIGSTRTH = zeros(OUTARRSIZE.l); 

LATGTH = zeros(OUTARRSIZE, 1 ); 

LATJTH = zeros(OUTARRSIZE, 1 ); 

DRUTH = zeros(OUT ARRSIZE, 1 ); 

DRURTH = zeros(OUTARRSIZE,l); 

DRUDEADTH = zeros(OUTARRSIZE,l); 

DRMPUUTH = zeros(OUTARRSIZE, 1); 

DECELTH = zeros(OUTARRSIZE.l); 

EXITPOSTH = zeros(OUTARRSIZE, 1); 

NEXITTH = zeros(OUTARRSIZE,l); 

VEXITTH = zeros(OUT ARRSIZE, 1 ); 

RUNTIMETH = zeros(OUTARRSIZE,l); 

LONGGTH = zeros(OUTARRSIZE, 1 ); 

LONGJTH = zeros(OUTARRSIZE, 1 ); 

POSEXITTH = zeros(OUTARRSIZE, 1 ); 

NEGEXITTH = zeros(OUTARRSIZE,l); 

ELEVTH = zeros(OUTARRSIZE, 1 ); 

DELRCTH = zeros(OUTARRSIZE, 1 ); 

DPERPTH = zeros(OUTARRSIZEJ); 

VPERPTH = zeros(OUTARRSIZE, 1 ); 

TIMETH = zeros(OUTARRSIZE, 1 ); 

FCTH = zeros(OUT ARRSIZE, 1 ); 

FNOSETH = zeros(OUTARRSIZE, 1 ); 

FMTH = zeros(OUTARRSIZE, I ); 

PMMTH = zeros(OUTARRSIZE,l); 

TRTH = zeros(OUTARRSIZE, 1 ); 

MNAEROTH = zeros(OUTARRSIZE,l); 

PSICTH = zeros(OUTARR.SIZE, 1 ); 

NCTH = zeros(OUTARRSIZE, 1 ); 

PSICGVTH = zeros(OUTARRSLZE, 1 ); 

THSTTH = zeros(OUTARRSIZE, 1 ); 

PBPERTH = zeros(OUTARRSIZE, 1 ); 

VWTH = zeros(OUTARRSIZE,l); 

VTWTH = zeros(OUTARRSIZE, 1 ); 

PSIPTH = zeros(OUTARRSIZE, 1 ); 

YAWTH = zeros(OUTARRSIZE, 1); 

MUNOSETH = zeros(OUTARRSIZE, 1 ); 

MUMAXTH = zeros(OUTARR.SIZE, 1 ); 
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MUMAXXTH = zeros(OUTARRSIZE, 1 ); 
MUMAXRTH = zeros(OUTARRSIZE,l); 
MUMAXLTH = zeros(OUTARRSIZE, 1 ); 
MUMAXCTH = zeros(OUTARRSIZE, 1 ); 
MUCTH = zeros(OUTARRSIZE, 1 ); 
MURTH = zeros(OUTARRSIZE, 1 ); 
MULTH = zeros(OUTARRSIZE, 1 ); 

FMRTH = zeros(OUTARRSIZE, 1 ); 

FMLTH = zeros(OUTARRSIZE, 1 ); 
BDRAGCTH = zeros(OUTARRSIZE,l); 
BDRAGRTH = zeros(OUTARRSIZE, 1 ); 
BDRAGLTH = zeros(OUTARRSIZE,l); 
PMCTH = zeros(OUTARRSIZE, 1 ); 

PMRTH = zeros(OUT A RRSIZE, 1 ); 

PMLTH = zeros(OUTARRSIZE,l); 
YNAVTH = zeros(OUTARRSIZE,l); 
XNAVTH = zeros(OUTARRSIZE,l); 
DECELCOMTH = zeros(OUTARRSIZE,l); 
DRAGCOMCTH = zeros(OUTARRSIZE, 1 ); 
RMPRATETH = zeros(OUTARRSIZE, 1 ); 
EABTH = zeros(OUTARRSIZE,l); 

PBTH = zeros(OUTARRSIZE, 1 ); 

SIGCTH = zeros(OUTARRSIZE, 1 ); 
YSTRIPETH = zeros(OUTARRSlZE, 1); 
ZEROTH = zeros(OUTARRSIZE,l); 

PBRTH = zeros(OUTARRSIZE, 1 ); 

PBLTH = zeros(OUTARRSIZE,l); 

FCRTH = zeros(OUTARRSIZE, 1 ); 

FCLTH = zeros(OUTARRSIZE,l); 
PSIMRTH = zeros(OUTARRSIZE, 1 ); 
PSIMLTH = zeros(OUTARRSIZE, 1 ); 
DRAGCOMLTH = zeros(OUTARRSIZE,l); 
DRAGCOMRTH = zeros(OUTARRSIZE,l); 
DRAGMAXLTH = zeros(OUTARRSIZE,l); 
DRAGMAXRTH = zeros(OUTARRSIZE,l); 
RADIUSTH = zeros(OUTARRSIZE, 1 ); 
TEMPI = zeros(OUTARRSIZE,l); 

TEMP2 = zeros(OUT ARRSIZE, I ); 

TEMP4 = zeros(OUT ARRSIZE, 1 ); 

TEMP5 = zeros(OUTARRSIZE, 1 ); 

TEMP6 = zeros(OUTARRSIZE, 1 ); 

TEMP3 = zeros(OUTARRSIZE, 1); 
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GLOB ALFILE. M is a file which declares as global so that they can be passed to a script 
file function in a common block rather than as a parameter and is called by RUNROTO 
once. 

global SP 
global SPM 
global WETMUT 
global WETMU 
global ICEMUT 
global 1CEMU 
global SNOWMUT 
global SNOWMU 
global SLUSHMUT 
global SLUSHMU 
global FLOODMUT 
global FLOODMU 
global CONSTMUT 
global CONSTMU 
global CONST2MUT 
global CONST2MU 

global PATCHNOW FIRSTPATCH LASTPATCH PATCHLEN PATCHTYPE PATCHOFTEN 
global DT “ 

global TMGD TNGD USEMUABORT 

global COMPUP DECLIM SST K2 RB OD HS DELS RP WS S MUROLL D1 K1 K2 DELB 
global spiralexitx spiralexity spiralexitx22 spiralexity22 spiralexitr22 
global TIMEOCCLAG FLATMU 

global tmpfyl 1 1 tmpfy222 mulimllst mulim21st MUSKIDT MUSKID 
global SW RWFC A B BC W 

global HCG ASEFF VEAS VTW XSTART REVERSE AUTOREV 
global CDFWD CDAFT CDFWD CGAFT CGFWD PCTMAC CGFWD 

global THS1T THS1 THSRT THSR THSTIT THSTI THSNT THSN THSRNT THSRN fyl 1 1 
fy222 3 

global psisk2 vskid2 murat2 mulim2 MD1 1 CG 

global psiskl vskidl muratl mulirol thstmulast BRKBUF DRATE DTIME WINGSPAN 
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FUNCPREDICTABORT.M is a MATLAB script function file called from 
RUNROTO.M once, containing the exit prediction algorithm. The parameters passed to 
it are variables predicted at touchdown. The following files are additional script file 
functions called by this algorithm. 

function abortearly = funcpredictabort(upre,utpre,xnavpre,vws$pre,vexit,exitpos,nexit) 

% these inputs are predicted values at touchdown 
% in the simulation time=0 at maingear touchdown 

global SW RWFC W ASEFF REVERSE AUTOREV 

global CDFWD CDAFT CGAFT CGFWD PCTMAC DECLIM TMGD TNGD 

global USEMUABORT TIMEOCCLAG 

global thstmulast BRKBUF DRATE DTIME MD1 1 CG 

% initialize variables 
thstmulast=0; 
abortearly =0; 
decelpreon=0; 

runtimepre= 1 /(upre/xnavpre); 

Ienfcms=250; 
disfcms=0; 
lenx =50; 
futt=0; 
fudt=0; 

tempw = utpre*1.689-upre; 

futut=utpre; 

deccnstu=upre; 

deccnstx=xnavpre; 

deccnstbd=exitpos-xnavpre; 

futu=upre; 

futut=(futu+tempw)/l .689; 

thstpre=functhrust(REVERSE,AUTOREV,futt,futut,xnavpre,exitpos,0.,futu i 0.5); 
CDRAG=CDFWD+((CDAFT-CDFWD)/(CGAFT-CGFWD))*(PCTMAC-CGFWD); 
adragmu = (utpre A 2)/295.37*CDRAG*SW; 
itemp= 1 +xnavpre/lenfcms; 

% how does weight and eg error affect exit prediction 
wpre=W+0; 
cgpre=CG-K); 
if(MDl 1) 

if(wpre>480000)wpre=480000;end; 
if(wpre<340000)wpre=340000;end; 
if(cgpre>.34)cgpre=.34;end; 
if(cgpre<. 12)cgpre=. 12;end; 
else 

if(wpre> 1 28000)wpre= 1 28000;end; 
if(wpre<82000)wpre=82000;end; 
if(cgpre>.34)cgpre=.34;end; 
if(cgpre<-0.008)cgpre=-.008;end; 
end; 

% CALCULATE NEW cgpre VARIABLES BASED ON FRACTION CG 
if(MDl 1) 

cginpre= 1311 .947+CG*295.779; 
apre=(cginpre-473 .43 7)/ 12.; 
bpre=( 1 442-cginpre)/ 12.; 
bcpre=( 1 472.62-cginpre)/12.; 

hcgpre=(209 .32-(sqrt(cginpre A 2+(2 1 ) A 2))*sin(atan(2 iycginpre)+0.00 1 93))/l 2.; 
else 
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cginpre=885.547+CG* 158.5 12; 
apre=(cginpre-97.998)/12.; 
bpre=(967. 1 -cginpre)/12.; 
bcpre=0; 

hcgpre=(83.029-(sqrt(cginpre A 2-K5. 1 ) A 2))*sin(atan(-5. l/cginpre)+0.0 1 8))/ 12.; 
end; 

temp=apre+bpre; 

[mumax,mumaxx]=funcfc(RWFC(itemp),upre,xnavpre,temp); 

% drag contributions due to crosswind 
dragcrs=abs( vwsspre)/5 7.3*4; 

% drag contributions 

otherd=(( + thstpre-adragniu)/wpre)*32.2+dragcrs; 

tempdlast=0; 

tempd= (vexit A 2- futu A 2)/(exitpos-xnavpre)/2.0 ; 
if(nexit>=3 | -tempd>7.0) 

DRATE=8.0; 

else 

DRATE=6.0; 

end; 

tempd=0; 

timeocclag=TIMEOCCLAG; 

timepreocclag=timeocclag; 

decelprerate=0; 

lastprexnav=xnavpre; 

lastpredeccalc=0; 

% advance aircraft location to next friction measurement (every lenx feet) 
for i = 200:lenx : 1 1450 
% start after touchdown point 
if(i > xnavpre & i < exitpos) 
futulast=futu; 

% delay decel cmd after touchdown 
if(-TMGD+futt < TNGD) 
tempd=0; 
end; 

if(-tempd>DECLIM) 

tempd=-9; 

end; 

% rate limit decel cmd 

if(abs(-tempd-tempdlast) > 1.33*fudt) 
if(-tempd-tempdlast > 0.) 

tempdset=tempdlast+ 1.33* fudt; 
else 

tempdset=tempdlast- 1 .33*fudt; 
end; 
else 

tempdset=-tempd; 

end; 

tempdlast=tempdset; 

% advance ground speed for travel distance using decel cmd, drag & thrust 
if(i-xnavpre > lenx) 

fhtu=sqrt(futu * futu+2 *(-tempdset+otherd)*(lenx)); 
else 

ftitu=sqrt(ftitu*ftitu+2*(-tempdset+otherd)*(i-xnavpre)); 

end; 

% calculate desired speed profile ground speed 
if(decelpreon) 
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temp=deccnstu-(i+BRKBUF-deccnstx)*(deccnstu-vexit)/deccnstbd; 

else 

temp=0; 

end; 

% limit ground speed at or above speed profile ground speed 
if(temp < vexit) 
temp=vexit; 
end; 

if(futu < temp) 
futu=temp; 
end; 

% advance time for travel distance using average of ground speed & last value 
if(i-xnavpre > lenx) 
ftidt=lenx/(ftitulast+futu)*2; 
else 

fudt=(i-xnavpre )/(futulast+futu)*2; 
end; 

futt=futt+fudt; 

% this is how you would normally calculate time except that we limit 
% futu above to const vel line 
% futt^futt+Cfutu-fiitulastyC-tempdset); 

% calculate airspeed, assumes winds are constant 
futut=(fiitu+tempw)/l .689; 

% calculate thrust using airspeed & whether braking is engaged 
temp=i; 
if(decelpreon) 

thstmu=functhrust(REVERSE,AUTOREV,fiitt,futut,temp,exitpos, 1 .,futu,0.5); 
else 

thstmu=functhrust(RE VERSE, AUTOREV,ftitt,futut,temp,exitpos,0.,ftitu,0.5); 
end; 

% calculate aircraft drag along runway using airspeed 
adragmu = (fiitut A 2)/295.37*CDRAG*SW; 

% calculate friction available from surface at main gear (mumaxx) 

% lenx stepsize may be smaller than friction measurement spacing 
if(disfcms<=0) 
disfcms=ienfcms; 
itemp= 1 + i/1 en fern s ; 
temp=i; 

temp2=apre+bpre ; 

[mumax,mumaxx]=ftmcfc(RWFC(itemp),futu,temp,temp2); 

end; 

disfcms=disfcms-lenx; 

% calculate required deceleration 

tempd = ((vexit ) A 2- futu A 2)/(exitpos-i)/2.0; 
deceltemppre=-tempd; 

% calculate AVAILABLE friction fraction available at main gear using 
% aircraft parameters (gear loading & avg anti-skid eff) 
muavailaiT=mumaxx*ASEFF*... 

(apre+hcgpre*(-adragmuywpreH-tempd/32.2))/(apre+(8*bpre+2*bcpre)/10.); 

% calculate NEEDED friction taking into account aircraft drag and thrust 
% *2 increases predict accuracy 

muneedarr=(-tempd/32.2-(-thstmu4-adragmu)/wpre)*2; 
otherd=((+thstmu-adragmu)/wpre)* 32 .2+dragcrs ; 
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% engage braking logic 

temp=( (exitpos-i )*(-tempd -lastpredeccalc)/... 

( i-lastprexnav)-tempd ); 
decelprerate = temp; 
lastprexnav= i; 
lastpredeccalc=-tempd; 
timepreocclast=timepreocclag; 
timepreocc=(exitpos-i)/((vexit+futu)/2.2)+futt+runtimepre; 

% TIME FOR DECEL AT 7.5 
tempt= (futu-vexit)/7.5; 

% DISTANCE OF DECEL AT 7.5 

tempdis=-(vexit A 2- futu A 2)/2.0/7.5; 

% TIME REQUIRED 

temp=BRKBUF/vexit+(exitpos-i-BRKBUF-tempdis)/futu+tempt; 
% TIME REMAINING 

dtemp=DTIME-futt-runtimepre; 

timepreocclag=timepreocc^timepreocclag-tiniepreocc)*exp(-0.5); 

% set decel command to available friction if braking is engaged 
if((decelpreon) |... 

((decelprerate>DRATE) |(timepreocc<DTIME)|(temp>dtemp)|... 
((timepreocclag>timepreocclast)&(futt>6)))) 
if(~decelpreon) 
deccnstu=futu; 
deccnstx=i; 
deccnstbd=exitpos-i; 
end; 

decelpreon=l; 

tempd=-muavailarr*32.2; 

else 

tempd=0; 

end; 

end; 

end; 

% exit for loop if aircraft is past ROTO exit 
% ABORT EXIT if NEEDED friction > AVAILABLE friction 
ifrUSEMUABORT) 
abortearly=muneedarr > muavailarr; 

% [muneedarr,muavailarr,futu] 
end; 

% ABORT EXIT if required deceleration > allowed deceleration 
if((~abortearly)&(deceltemppre>DECLIM)) 
abortearly= 1 ; 
end; 
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FUNCTHRUST.M is a script function calculating forward and reverse thrust and is 
called by FUNCPREDICTABORT. 

function thst = functhrust(reverse,autorev,time,ut,xnav,exitpos,aurevmult,u,dttau) 

global THS1T THS1 THSRT THSR THSTIT THSTI THSNT THSN THSRNT THSRN 
thstmulast 

% auto-reverse thrust control logic 
% do not use auto-reverse thrust until brakes start 
if((~reverse)|((xnav+250)>exitpos)|(u<70* 1 .689)) 
if (time<8.0) 

thst=funcinterp(ut,THS 1 T,THS 1 , length(THS 1 T)); 
else 

thst=funcinterp(ut,THSNT,THSN,length(THSNT)); 

end; 

else 

% reverse thrust for non-dry runway 
if (time<THSTIT(l)) 

thst=fiincinterp(time ,THSTIT,THSTI,length(THSTIT)); 

% thrust arrays have to be set right for this to worldscale spool up 
if(ut<THSRT(2)) 

temp^funcinterp^t , THSRT, THSR,length(THSRT)); 
temp=temp/THSR(2); 
else 

temp= 1 ; 
end; 

thst=thst*temp; 

else 

thst=funcinterp(ut , THSRT, THSR,length(THSRT)); 
end; 

thsttemp=funcinterp(ut,THSRNT,THSRN,length(THSRNT)); 

% add reverse idle limit here 

if((th$t<thsttemp) & autorev & (time>5.)) 
temp= aurevmult*(thst-thsttemp)+thsttemp; 

% engine lag, first 1 0 sec already lagged 

thst = temp+(thstmulast -temp)*exp(-dttau ); 
end; 
end; 

thstmulast=thst; 
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FUNCFC.M is a script function which calculates surface friction and is called by 
FUNCPREDICTABORT. 

function [mumax,mumaxx]=funcfc(wet,U,x,nlen) 

global SP 
global SPM 
global WETMUT 
global WETMU 
global ICEMUT 
global ICEMU 
global SNOWMUT 
global SNOWMU 
global SLUSHMUT 
global SLUSHMU 
global FLOODMUT 
global FLOODMU 
global CONSTMUT 
global CONSTMU 
global CONST2MUT 
global CONST2MU 

global PATCHNOW FIRSTPATCH LASTPATCH PATCHLEN PATCHTYPE PATCHOFTEN 

munose= 1 ; 
mumain=l; 

% get nose gear mu 
% assume PS1P=0 

if(PATCHNOW&((x+nlen)>FIRSTPATCH)&((x+nlen)<LASTPATCH)& ... 

( fix( fix((x+nlen)/PATCHLEN)/PATCHOFTEN)= ... 

fix((x+nlen)/PATCHLEN)/PATCHOFTEN)) 

wetpatch=PATCHTYPE; 

else 

wetpatch=wet; 

end; 

% dry runway 

mumax = (0.93 - 0.001 1* SP )*(!. - 0.0013*( U/1.689)); 

% wet runway 

if(wetpatch=l) 

munose=funcinterp(U/ 1.689, WETMUT, WETMU, length(WETMUT)); 

% icy runway 
else 

if(wetpatch=2) 

munose=funcinterp(U/1.689,ICEMUT,ICEMU,length(ICEMUT)); 

% snowy runway 
else 

if(wetpatch=3) 

munose=funcinterp(U/l. 689, SNOWMUT, SNOWMU,length(SNOWMUT)); 

% slushy runway 
else 

if(wetpatch==4) 

munose=funcinterp(U/1.689,SLUSHMUT,SLUSHMU,length(SLUSHMUT)); 

% flooded runway 
else 

if(wetpatch=5) 

munose=funcinterp(U/ 1 .689,FLOODMUT,FLOODMU,length(FLOODMUT)); 

% const runway , read friction coefficient in from input directly 
else 

if(wetpatch=6) 

munose=funcinterp(U/l. 689, CONSTMUT, CONSTMU, length(CONSTMUT)); 
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else 

if(wetpatch==7) 

munose=funcinterp(U/1 .689,CONST2MUT,CONST2MU,length(CONST2MUT)); 
end; 
end; 
end; 
end; 
end; 
end; 
end; 

% get main gear mu 

if(PATCHNOW&((x)>FIRSTPATCH)&((x)<LASTPATCH)& ... 

( fix( fix((x)/PATCHLEN)/PATCHOFTEN)= ... 

fix((x)/PATCHLEN)/PATCHOFTEN)) 

wetpatch=PATCHTYPE; 

else 

wetpatch=wet; 

end; 

% dry runway 

mumaxx= (0.93 - 0.001 1* SPM)*(1. - 0.0013*( U/1.689)); 
if(wetpatch=l) 

mumain=funcinterp(U/ 1 .689, WETMUT, WETMUJength(WETMUT)); 

% icy runway 
else 

if(wetpatch— 2) 

mumain=funcinterp(U/1.689,ICEMUT,ICEMU,length(ICEMUT)); 

% snowy runway 
else 

if(wetpatch=3) 

mumain=funcinterp(U/ 1 .689,SNOWMUT,SNOWMU,length(SNOWMUT)); 

% slushy runway 
else 

if(wetpatch==4) 

mumain=funcinterp(U/1.689,SLUSHMUT,SLUSHMU,length(SLUSHMUT)); 

% flooded runway 
else 

if(wetpatch=5) 

mumain=funcinterp(U/1.689,FLOODMUT,FLOODMU,length(FLOODMUT)); 

% const runway , read friction coefficient in from input directly 
else 

if(wetpatch==6) 

mumain=funcinterp(U/1.689,CONSTMUT,CONSTMU,length(CONSTMUT)); 

else 

if(wetpatch==7) 

mumain=funcinterp(U/l,689,CONST2MUT,CONST2MU,length(CONST2MUT)); 

end; 

end; 

end; 

end; 

end; 

end; 

end; 

if(wetpatch>5) 

mumaxx=mumain; 

mumax=munose; 

else 

mumaxx=mumaxx*mumain; 

mumax=mumax*munose; 


A-33 



end; 
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FUNCINTERPY.M is a script function returning the y value of the centerline and is called 
by FUNCPREDICTABORT. 

function y=funcinterpy(x,exitpos) 

global spiralexitx spiralexity 

newx=x-exitpos; 
if (newx<0) 

y-o.; 

else 

if (newx>2200.688) 
y=600.; 
else 

if (newx< 1362.1) 

% dispC Point 1, funcinterpy'),keyboard 

y=funcinterp(newx,spiralexitx,spiralexity,length(spiralexitx)); 

else 

if (newx<l 800.5 17) 

y=-460.397+0.5*(newx-149.669)/0. 8660254; 
else 

endcrv=2200; 

smrad=800; 

a=l.; 

b* 2. *200.; 

c=newx*newx-2.*endcrv*newx+endcrv*endcrv-smrad*smrad+200.*200.; 

y=(-b+sqrt(b*b-4.*a*c))/(2.*a); 

end; 

end; 

end; 
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FUNCINTERP.M is a script function interpolating vectors and is called by 
FUNCPREDICTABORT. 

function y = funcinterp(x,xarray,yarray,npts) 
if (xarray(l) < xarray(npts)) 
for i=2:npts 

if (xarray(i) > x)|(i=npts) 

y = (x-xarray(i-l))/(xarray(i)-xaiTay(i-l))*(yarray(i)-yaiTay(i-l)) + yarray(i-l); 

break; 

end; 

end; 

else 

for i=npts-l:-l:I 
if(xarray(i) > x)|(i=l) 

y=(x-xarray(i+ 1 ))/(xarray(i)-xarray(i+ 1 )) # (yarray(i)-yarray(i+ 1 ))+yarray(i+ 1 ); 

break; 
end; 
end; 
end ; 
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FUNCPERPDIST.M is a script function returning the aircraft’s perpendicular distance to 
the centerline and is called by FUNCPREDICTABORT. The XPATH input is already 
determined to be the X location of the centerline perpendicular to the aircraft, 
function dist = funcperpdist(x,y,xpath,exitpos) 

% REAL*8 DIST 


ypath = funcinterpy(xpath,exitpos); 
if(abs(y-ypath)~=0) 

dist=sqrt((x-xpath) A 2+(y-ypath) A 2)*abs(y-ypath)/(y-ypath); 

else 

dist=abs(x-xpath); 

end; 
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SCRIPTWRITEPLOTDATA.M is a script file to create and write output variables from 
the simulation run and is called by RUNROTO. The following files are used by 
scriptwriteplotdata. 

global fid 

i=TIME/INTVL+ 1 ; 
if(i> 1 000)i= 1000;end; 
for TSTEP=l:i 

skip=((TSTEP- 1 )/DT* INTVL)+ 1 ; 

% create variables that are created from others before 
% creating others 


% DISPLAY WHEN WING TIP EXITS RUNWAY SIDE, VARIABLE ON RUNWAY WIDTH=75 
clrdist=WINGSPAN*cos(PSITH(skip)/57.3)+75; 
if(abs(YTH(skip))>clrdist) 

RUNT1METH(TSTEP)=RUNTIMETH(TSTEP -1); 
else 

RUNTIMETH(TSTEP)=RUNTIMETH(skip); 

end; 

RUNTIMETH( 1 )= 1 /((VEAS+VTW)* 1 .689/XSTART); 
if(LONGJTH(TSTEP)>0. 1 5)LONGJTH(TSTEP>=0. 1 5;end; 
if(LONGJTH(TSTEP)< -0.45)LONGJTH(TSTEP)= -0.45;end- 


EXITLOC( 1 ),EXITPOSTH(skip)); 
else 


if((XTH(skip)<800+EXITLOC(2))&(XTH(skip)>EXITLOC(2))&(XTH(skip)<EXITPOSTH(skip))) 

EXITLO^xTl^^^ 

else 

, ( ( XTH(!ldp)^OWKirilOW))*0™^^EmWC^^n^^^XITPOSTH(!kip))) 

EXITLOC(3),EXITPOSTH(skip»; 1 P 

else 

YSTRIPETH(TSTEP)=funcinterpy(XTH(skip),EXITPOSTH(skip)); 

end; 

end; 

end; 

TEMF=XTH(skip)-EXITPOSTH(skip); 

if(TEMP<0.0) 

POSEXITTH(TSTEP) = (( 1 50-TREAD)/2. -LATBIAS); 
else 

if (TEMP<680 ) 

POSEXITTH(TSTEP) = (((90+60*(680 -TEMP)/680 ) -TREADV2. -LATBIAS)- 
else 

POSEXITTH(TSTEP) = ((90-TREAD)/2. -LATBIAS); 
end; 
end; 

NEGEXITTH(TSTEP) = -POSEXITTH(TSTEP) ; 

ZEROTH(TSTEP) =0; 
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XTH(TSTEP) = XTH(skip); 

PB PERTH(TSTEP) =(PBTH(skip) +PBRTH(skip) +PBLTH(skip) )/3./PSUP*100.; 
MUMAXRTH(TSTEP) =MUMAXXTH(skip) ; 

MUMAXLTH(TSTEP) = -MUMAXXTH(skip) ; 

MUMAXCTH(TSTEP) =MUMAXXTH(skip) ; 

MUCTH(TSTEP) =sqrt(FCTH(skip) A 2+BDRAGCTH(skip) A 2)/PMCTH(skip) ; 
MURTH(TSTEP) =sqrt(FMRTH(skip) A 2+BDRAGRTH(skip) A 2)/PMRTH(skip) ; 
MULTH(TSTEP) = -sqrt(FMLTH(skip) A 2+BDRAGLTH(sklp) A 2)/PMLTH(skip) ; 

TIMETH(TSTEP)=TIMETH(skip); 

DXTH(TSTEP) = DXTH(skip); 

DUTH(TSTEP) = DUTH(skip); 

UTH(TSTEP) = UTH(skip)/ 1.689; 

YTH(TSTEP) = YTH(skip); 

DYTH(TSTEP) = DYTH(skip); 

DVTH(TSTEP) = DVTH(skip); 

PSITH(TSTEP) = PSITH(skip); 

DPSITH(TSTEP) = DPSITH(skip); 

DDPS ITH(TSTEP) = DDPSITH(skip); 

SIGTH(TSTEP) = SIGTH(skip) ; 

SIGSTRTH(TSTEP) = SIGSTRTH(skip); 

LAT GTH(TSTEP) = LATGTH(skip); 

LATJTH(TSTEP) = LATJTH(skip) ; 

DRUTH(TSTEP) = DRUTH(skip); 

DRURTH(TSTEP) = DRURTH(skip) ; 

DRUDEADTH(TSTEP) = DRUDEADTH(skip) ; 

DRMPUUTH(TSTEP) = DRMPUUTH(skip); 

DECELTH(TSTEP) = DECELTH(skip); 

NEXITTH(TSTEP) = NEXITTH(skip); 

VEXITTH(TSTEP) = VEXITTH(skip); 

RUNTIMETH(TSTEP) = RUNTIMETH(skip); 

LONGGTH(TSTEP) = LONGGTH(skip); 

LON G JTH(TSTEP) = LONGJTH(skip); 

ELEVTH(TSTEP) = ELEVTH(skip) ; 

DELRCTH(TSTEP) = DELRCTH(skip) ; 

DPERPTH(TSTEP) = DPERPTH(skip); 

VPERPTH(TSTEP) = VPERPTH(skip); 

FCTH(TSTEP) = FCTH(skip); 

FNOSETH(TSTEP) = FNOSETH(skip); 

FMTH(TSTEP) = FMTH(skip) ; 

PMMTH(TSTEP) = PMMTH(skip); 

TRTH(TSTEP) = TRTH(skip); 

MNAEROTH(TSTEP) = MNAEROTH(skip) ; 

PSICTH(TSTEP) = PSICTH(skip); 

NCTH(TSTEP) = NC; 

PSICGVTH(TSTEP) = PSICGVTH(skip) ; 

THSTTH(TSTEP) = THSTTH(skip); 

VWTH(TSTEP) = VWTH(skip); 

VTWTH(TSTEP) = VTWTH(skip); 

PSIPTH(TSTEP) = PSIPTH(skip); 

YAWTH(TSTEP) = YAWTH(skip); 

MUNOSETH(TSTEP) = MUNOSETH(skip); 

MUMAXTH(TSTEP) = MUMAXTH(skip) ; 

FMRTH(TSTEP) = FMRTH(skip); 

FMLTH(TSTEP) = FMLTH(skip); 

BDRAGCTH(TSTEP) = BDRAGCTH(skip); 

BDRAGRTH(TSTEP) = BDRAGRTH(skip); 

BDRAGLTH(TSTEP) = BDRAGLTH(skip); 
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PMCTH(TSTEP) = PMCTH(skip); 

PMRTH(TSTEP) = PMRTH(skip); 

PMLTH(TSTEP) = PMLTH(skip); 

YNAVTH(TSTEP)=YNAVTH(skip); 

XN A VTH(TSTEP) - XNAVTH(skip); 

DECELCOMTH(TSTEP) = DECELCOMTH(skip); 

DRAGCOMCTH(TSTEP) = DRAGCOMCTH(skip); 

RMPRATETH(TSTEP) = RMPRATETH(skip); 

EABTH(TSTEP) = EABTH(skip); 

PBTH(TSTEP) = PBTH(skip); 

SIGCTH(TSTEP) = SlGCTH(skip); 

PBRTH(TSTEP) = PBRTH(skip) ; 

PBLTH(TSTEP) = PBLTH(skip); 

RADIUSTH(TSTEP) = RADIU STH(skip); 
end; 

TSTEP=i; 

% PRINT OUT TIME HISTORY ARRAYS AND CROSS PLOT ARRAYS 
fprint^fid,' ********* OUTPUT DATA ************** \n')- 

fprintf(fid,' X Y \n'); 

funccrossplot(XTH,YTH,TSTEP); 
fprintf(fid,' X YSTRIPE \n'); 

funccrossplot(XTH, Y STRIPETH,TSTEP); 
fprintf(fid,' TIME U Y PSI YSTRIPE \n'); 

functimhist(TIMETH,UTH,YTH,PSITH,YSTRIPETH,TSTEP); 
fprintf(fid,' X DV U PSI YSTRIPE \n'); 

functimhist(XTH,DVTH,UTH,PSITH,YSTRIPETH,TSTEP); 
fprintflfid,’ TIME DV DY DPSI DU \tf); 

functimhist(TIMETH,DVTH,DYTH,DPSITH,DUTH,TSTEP); 
fprintf(fid,' TIME DV SIGSTR DELRC DPERP \n'); 

functimhist(TIMETH,DVTH,SIGSTRTH,DELRCTH,DPERFTH,TSTEP)- 
fprintf(fid,' TIME DDPSI FC FNOSE FM W); 

functimhist(TIMETH,DDPSITH,FCTH,FNOSETH,FMTH,TSTEP); 
fprintf(fid,' TIME PMM PSIP THST TR \n'); 

functimhist(TIMETH,PMMTH,PSIPTH,THSTTH,TRTH,TSTEP); 
fprintf(fid,' TIME MNAERO PSIC NC PSICGV \n’); 

functimhist(TIMETH,MNAEROTH,PSICTH,NCTH,PSICGVTH,TSTEP); 
fprintf(fid,' X MUNOSE MUMAX DPERP VPERP \n’); 

functimhist(XTH,MUNOSETH,MUMAXTH,DPERPTH,VPERPTH,TSTEP); 
fprintf(fid,' X LATG SIG YAW DV \n'); 

ftjnctimhist(XTH,LATGTH,SIGTH,YAWTH,DVTH,TSTEP); 
fprintf^fid,’ SIGSTR FNOSE LATG MUNOSE DV\n*); 
functimhist(SIGSTRTH,FNOSETH,LATGTH,MUNOSETH,DVTH,TSTEP); 
fprintf(fid,' X DRU DRUR DRUDEAD DRMPUU \n’); 
functimhist(XTH,DRUTH,DRURTH > DRUDEADTH,DRMPUUTH,TSTEP); 
lprintf(fid,’ X DECELCOM NEXIT RUNWAYTIME LONGG \n'); 
functimhist(XTH,DECELCOMTH,NEXITTH,RUNTIMETH,LONGGTH,TSTEP); 
fprintf(fid,' X POSEXIT NEGEXIT LATJ LONGJ \n'); 
functimhist(XTH,POSEXITTH,NEGEXITTH > LATJTH,LONGJTH,TSTEP); 
fprintf(fid,' XPERCENT BRK THRUST TAILWIND CROSSWIND\n'); 
functimhist(XTH,PBPERTH,THSTTH,VTWTH,VWTH,TSTEP); 
fprintf(fid,' X XNAV YNAV ELEV CROSSWIND\n*); 
functimhist(XTH,XNAVTH,YNAVTH,ELEVTH,VWTH,TSTEP); 
fprintf(fid,' X MUMAXR MUR MUMAXC MUC \n'); 

functimhist(XTH,MUMAXRTH,MURTH,MUMAXCTH,MUCTH,TSTEP); 
fprintf(fid,' X MUMAXL MUL DX ZERO \n0; 

flinctimhist(XTH,MUMAXLTH,MULTH,DXTH»ZEROTH,TSTEP); 

fclose(fId); 
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FUNCCROSSPLOT.M is a script function write format and is called by 

SCRIPTWRITEPLOTDATA. 

function funccrossplot(paral,para2,tstep) 

global fid 


fprintf(fid,'PLOTDEF\n'); 
fprintf(fid,’XREAD( 1 )= 1 \n’); 
fprintf(fid.’YREAD( 1 )=2\n'); 
fprintf(fid,'END\n’); 
fprintf(fid,'DATA\n'); 
for i=l:tstep 

fprintf(fid,'% 1 3 .6e % 1 3 ,6e\n',para 1 (i),para2(i)); 
end; 

fprintf(fid,'END\n'); 
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FUNCTIMEHIST.M is a script function write format and is called by 
SCRIPTWRITEPLOTDATA. 

function functimhist(tim,paral ,para2,para3,para4,tstep) 


global fid 

fprintf(fid,'PLOTDEF\n'); 
fprintf(fld,'XREAD( 1 )=1 \n'); 
fprintflfid,' YRE AD( 1 )=2\n’); 
fprintf(fid,'YREAD(2)=3\n’); 
fprintf(fid,'YREAD(3)=4\n'); 
fprintf(fld,'YREAE>(4)=5\n'); 
fprintf(fid,'END\n’); 
f^rintf(fid.'DATAn'); 
for i=l: tstep 

fprintf(fid,'% 1 3 .6e %13.6e %13.6e %13.6e %13.6e\n , ,tim(i),paral(i),para2(i),para3(i),para4(i)); 
end ; 

fprintf(fid,’END\n'); 
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FUNCINPUT.C is C function file which is converted to a CMEX file and called by 
ROTO SIMULINK picture AERO in file MODELROTOBIG.M. It returns the aircraft 
lateral deviation from the centerline and aircraft heading. It uses CENTERLINE.M for 
spiral exit centerline coordinates (2 feet apart). 

#include <math.h> 

#include <stdio.h> 

#include "mex.h" 

/* Input Arguments */ 

#define UP_IN prhs[0] 

/* Output Arguments */ 

#define YP_OUT plhsfO] 

#defme inputvec 4 


double funcinterp(double x,double xarray[],double yarray[],int npts) 

{ 

double y; 
int i; 


y=-1.0; 

if (xarrayfO] < xarray[npts-l]) 

{ 

for (i=l;i<“npts-l;i++) 

{ 

if ((xarray[i] > x) ||(i==npts-l)) 

{ 

y = (x-xarray[i- 1 ])/(xarray [i]-xairay[i- 1 ])*(yarray [i]-yarray[i- 1 ]) + yarray[i-l]; 
return y; 

}; 

}; 

} 

else 

{ 

for (i=npts-2;i>=0;i — ) 

{ 

if((xarray[i] > x)||(i==0)) 

{ 

y=(x-xarray[i + l])/(xarray[i]-xarray[i+l])*(yarrayti]-yarray[i+l])-t-yarray[i+l]; 
return y; 

}; 

}; 

} 

return y; 

} 


double funcinterpy(double x, double exitpos) 

{ 

double y; 
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double newx,endcrv,smrad,a,b,c; 

double *spiralexitx,*spiralexity; 
int lenspiralexit; 

spiralexitx=mxGetPr(mexGetGlobal("spiralexitx")); 

spiralexity=mxGetPr(mexGetGlobal("spiralexity")); 

lenspiralexit=mxGetM(mexGetGlobaI("spiralexitx")); 


newx=x-exitpos; 
if (newx<0) 

{ 

y=0.; 

} 

else 

{ 

if (newx>2200.688) 

{ 

y=600.; 

} 

else 

{ 

if (newx<l 362.1) 

{ 

y=funcinterp(newx,spiralexitx,spiralexity, lenspiralexit); 

} 

else 

{ 

if (newx<l 800.5 17) 

{ 

y=-460.397+0.5*(newx- 1 49.669)/0.8660254; 

} 

else 

{ 

endcrv=2200; 

smrad=800; 

a=l.; 

b= 2. *200.; 

c=newx*newx-2.*endcrv*newx+endcrv*endcrv-smrad*smrad+200.*200.; 

y=(-b+sqrt(b*b-4.*a*c))/(2.*a); 

}; 

}; 

}; 

} 

return y; 

} 


double funcperpdist(double x.double y, double xpath,double exitpos) 

{ 

double dist,ypath; 
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ypath = funcinterpy(xpath,exitpos); 
if(fabs(y-ypath) ! =0) 

{ 

dist=sqrt(pow((x-xpath),2}+pow((y-ypath),2))*fabs(y-ypath)/(y-ypath); 

} 

else 

{ 

dist=fabs(x-xpath); 

} 

return dist; 

} 


/* function [dperp,vperp,psip] = funcinput(u) */ 
void funcinputfdoubie yout[], double u[]) 

{ 

double x,y,exitpos,psi; 
double olddp; 
double DT; 

double fpathx,dlat,far,near,mid,fardst,neardstjniddst,dperp,vperp; 
double lastmiddst,xclose,yclose,xplus,yplus,psip,dyp,dxp; 


DT=*mxGetPr(mexGetGlobal("DT")); 

/* olddp=*mxGetPr(mexGetGlobal(”olddp"));*/ 

x=u[0]; 

y=u[l]; 

exitpos=u[2]; 

psi=u[3]; 

/* printf(”%f %f %f %f\n",x,y,exitpos,psi);*/ 

/* AIRCRAFT DISPLACEMENT PERPENDICULAR TO PATH, CLOSEST POINT ON PATH */ 
fpathx= funcinterpy(x,exitpos); 
dlat=fabs(y-fpathx); 
if(dlat< 1.0) 

{ 

dlat=1.0; 

}; 

far=(x+dlat); 

near=(x-dlat); 

mid=x; 

fardst=funcperpdist(x,y ,far,exitpos) ; 
fardst=fabs(fardst); 

neardst=funcperpdist(x,y,near,exitpos); 
neardst=fabs(neardst); 
middst=funcperpdist(x,y,mid,exitpos); 
dperp= 10000; 
lastmiddst=0.0; 

while (((fabs(fabs(middst)-fabs(dperp))) > 0.0000005) && (middst != lastmiddst)) 

{ 

dperp=middst; 
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if (fardst-neardst < 0) 

{ 

near=mid; 

neardst=fabs(middst); 

} 

else 

{ 

far=mid; 

fardst=fabs(middst); 

} ; 

mid=(far+near)/2.0; 

lastmiddst=middst; 

middst = funcperpdist(x,y,mid,exitpos); 

}; /*% end of while loop*/ 

dperp=middst; 

xclose=mid; 

yclose=funcinterpy(xclose,exitpos); 

/*% AIRCRAFT VELOCITY PERPENDICULAR TO PATH 
cannot take derivatives in cmex file, file called more than once per iteration. 
vperp=(dperp-o!ddp)/DT ; 
olddp=dperp; 

mxSetPr(mexGetGlobaI("olddp"),&olddp);*/ 

/*% AIRCRAFT YAW ANGLE WITH RESPECT TO PATH*/ 
dxp=.00 1 ; 
xplus=xclose+dxp; 
yplus=funcinterpy(xplus,exitpos); 
dyp=yplus-yclose; 
psip=psi-57.3 *atan(dyp/dxp); 


yout[0]=dperp; 

/* yout[l]=vperp;*/ 
yout[I]=psip; 


} 


void mexFunction( 


int 


nlhs. 

Matrix 

*plhs[], 


int 


nrhs. 

Matrix 

) 

*prhs[] 


double 

*yp; 


double 

*u; 


unsigned int 

m,n; 


/* Check for proper number of arguments */ 
if(nrhs !=1) { 

mexErrMsgTxt(”Funcinput routine requires one vector input argument.")- 
} else if (nlhs > 1) { 

mexErrMsgTxt("Funcinput routine produces one vector output argument."); 
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/* Check the dimensions of u. */ 

m = mxGetM(UPIN); 
n = mxGetN(UPIN); 

if (!mxIsNumeric(UP_rN) || mxIsComplex(UP_IN) || 
!mxIsFull(UP_IN) || !mxIsDouble(UP_IN) || 

(!(m = inputvec && n — 1) &&!(m = 1 && n = inputvec)) ) 

mexErrMsgTxt("Funcinput routine requires 4 input vector."); 

} 

/* Create a Matrix for the return argument */ 

YP OUT = mxCreateFull(l,2, REAL); 

/* Assign pointers to the various parameters */ 

yp = mxGetPr(YP_OUT); 
u = mxGetPr(UP_IN); 

/* Do the actual computations in a subroutine */ 

funcinput(yp,u); 

return; 

} 
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FUNCINPUTNAV.C is C function file which is converted to a CMEX file and called by 
ROTO SIMULINK picture CTRLLAWS20HZ in file MODELROTOBIG.M. It returns 
the aircraft lateral deviation from the centerline and aircraft heading. It uses 
CENTERLINE22.M for spiral exit centerline coordinates. (~100 feet apart) 

#include <math.h> 

#include <stdio.h> 

#include "mex.h" 

/* Input Arguments */ 

#define UP_IN prhs[0] 

/* Output Arguments */ 

#define YP_OUT plhs[0] 

#define inputvec 4 


double funcinterpnav(double x, double xarray[], double yarray[],int npts) 

double y; 
int i; 


y=-1.0; 

if (xanay[0] < xarray[npts~l]) 

{ 

for (i=l;i<=npts-l;i++) 

{ 

if ((xarrayfi] > x)||(i=npts-l)) 

{ 

y = (x-xarray[i-l])/(xairay[i]-xarray[i-l])*(yarrayti]-yarray[i-lj) + yanay[i-l]; 

return y; 

}; 

}; 

} 

else 

{ 

for (i=npts-2;i>=0;i — ) 

{ 

if((xarray[i] > x) ||(i==0» 

{ 

y=(x-xarray [i+ 1 ])/(xarray[i]-xarray[i+ 1 ])*(yarray[i]-yarray[i+ 1 ])+yanay[i+ 1 ]; 
return y; 

}; 

}; 

} 

return y; 

} 


double funcintrpynav(double x. double exitpos) 

{ 
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double y; 

double newx,endcrv ? smrad,a t b,c; 

double* spiralexitx22,*spiralexity22; 
int lenspiralexit22; 

spiralexitx22=nlxGetP^(mexGetGlobaI("spi^alexitx22 ,, )); 

spiralexity22=mxGetPr(mexGetGlobal( ,, spiralexity22 M )); 

lenspiralexit22=mxGetM(mexGetGlobal( M spi^alexitx22 , ')); 


newx=x-exitpos; 
if (newx<0) 

{ 

y=0.; 

} 

else 

{ 

if (newx>2200.688) 

{ 

y=600.; 

} 

else 

{ 

y=funcinterpnav(newx,spiralexitx22,spiralexity22,lenspiralexit22); 

}; 

}; 


return y; 

} 


double funcperpdistnav(double x,double y,double xpath, double exitpos) 


double dist,ypath; 


ypath = funcintrpynav(xpath, exitpos); 
if(fabs(y-ypath)!=0) 

{ 

dist=sqrt(pow((x-xpath),2)+pow((y-ypath),2))*fabs(y-ypath)/(y-ypath); 

} 

else 

{ 

dist=fabs(x-xpath); 

}; 

return dist; 

} 
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!* function [dperp,vperp,psip] = funcinputnav(u) */ 
void funcinputnav(double yout[], double u[]) 

{ 

double x,y,cgx,cgy,exitpos,psi; 
double olddpl,A,B; 
double DT; 

double fpathx,diat,far,near,mid,fardst,neardst,middst,dperp,vperp; 
double lastmiddst,xclose,yclose,xplus,yplus,psip,dyp,dxp; 

DT=*mxGetPr(mexGetGlobaI("DT")); 

/* olddpl =*mxGetPr(mexGetGlobal("olddpl "));*/ 
A=*mxGetPr(mexGetGlobal("A")); 
B=*mxGetPr(mexGetGlobal("B")); 

cgx=u[0]; 

cgy=u[l]; 

exitpos=u[2]; 

psi=u[3]; 

x=cgx+(A+B)*cos(psi/57.3); 
y=cgy+(A+B)*sin(psi/57.3 ); 


/* AIRCRAFT DISPLACEMENT PERPENDICULAR TO PATH, CLOSEST POINT ON PATH */ 
fpathx= funcintrpynav(x,exitpos); 
dlat=fabs(y-fpathx); 
if(dlat < 1 .0) 

{ 

dlat=1.0; 

}; 

far=(x+dlat); 

near=(x-dlat); 

mid=x; 

fardst=funcperpdistnav(x,y,far,exitpos); 

fardst=fabs(fardst); 

neardst=funcperpdistnav(x,y,near,exitpos); 

neardst=fabs(neardst); 

middst=funcperpdistnav(x,y,mid,exitpos); 

dperp= 1 0000; 

lastmiddst=0.0; 

while (((fabs(fabs(middst)-fabs(dperp))) > 0.0000005) && (middst != lastmiddst)) 

dperp=middst; 
if (fardst-neardst < 0) 

{ 

near=mid; 

neardst=fabs(middst); 

} 

else 

{ 

far=mid; 

fardst=fabs(m iddst); 

} ; 

mid=(far+near)/2.0; 

lastmiddst=middst; 

middst = funcperpdistnav(x,y,mid,exitpos); 

}; !*% end of while loop*/ 
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dperp=middst; 

xclose=mid; 

yclose=funcintrpynav(xclose,exitpos); 

/*% AIRCRAFT VELOCITY PERPENDICULAR TO PATH 
DT would be COMPUP if it were running at 20 hz 

cannot create derivatives in mex files because it is called more than once per iteration 
vperp=(dperp-olddp 1 )/DT ; 
olddpl=dperp; 

mxSetPr(mexGetGlobal( M olddpr , ),&olddpl);*/ 

/*% AIRCRAFT YAW ANGLE WITH RESPECT TO PATH*/ 
dxp=.001 ; 
xpius=xclose+dxp; 
yplus=funcintrpynav(xplus,exitpos); 
dyp=yplus-yclose; 
psip=psi-57.3*atan(dyp/dxp); 


yout[0]=dperp; 

/* youtflj^vperp;*/ 
yout[l]=xclose; 

} 


void mexFunction( 


int 


nlhs, 

Matrix 

*plhs[], 


int 


nrhs, 

Matrix 

) 

*pihs[] 


double 

*yp; 


double 

*u; 



unsigned int m,n; 

/* Check for proper number of arguments */ 
if(nrhs !=1) { 

mexErrMsgTxt("Funcinputnav routine requires one vector input argument."); 

} else if (nlhs > 1) { 

mexErrMsgTxt("Funcinputnav routine produces one vector output argument.”); 

} 


/* Check the dimensions of u. */ 

m = mxGetM(UPIN); 
n = mxGetN(UP _IN); 

if (!mxlsNumeric(UP_IN) || mxIsComplex(UP_IN) || 
!mxlsFull(UP_IN) || !mxIsDouble(UP_IN) || 

(!(m = inputvec && n == 1) &&!(m = 1 && n = inputvec)) ) 

mexErrMsgTxt("Funcinputnav routine requires 4 input vector.”); 

} 

/* Create a Matrix for the return argument */ 
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YPOUT = mxCreateFull(l,2, REAL); 

/* Assign pointers to the various parameters */ 

yp = mxGetPr(YP_OUT); 
u = mxGetPr(UP_JN); 

/* Do the actual computations in a subroutine */ 

funcinputnav(yp,u); 

return; 
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FUNCINTERPR.C is C function file which is converted to a CMEX file and called by 
ROTO SIMULINK picture CTRLLAWS20HZ in file MODELROTOBIG.M. It returns 
the centerline radius and uses CENTERLINE22.M. (-100 feet apart) 


#include <math.h> 

#include <stdio.h> 

#include "mex.h" 

/* Input Arguments */ 

#define UPIN prhs[0] 

/* Output Arguments */ 

#define YP_OUT plhs[0] 

#define inputvec 3 


double interpr(double x,double xarrayf], double yarray[],int npts) 

{ 

double y; 
int i; 


y=-10; 

if (xarray[0] > xarray[npts-l]) 

{ 

for (i= 1 ;i<=npts- 1 ;i++) 

{ 

if ((xarray[i] < x)||(i==npts- 1 )) 

{ 

/* This logrithmic interpolation technique tried to represent the spiral nature of the exit radius. 

Without it, the interpolated radius would provide more lead.*/ 

y = -(pow( 1 0, 1 -((x-xarray [i- 1 ])/(xarray[i]-xarray[i- 1 ])))- 1 .0)/9.0*(yarray[i]-yarray[i- 1 ]) + yarray [i]; 
return y; 

}; 

}; 

} 

else 

{ 

for (i=npts-2;i>=0;i— ) 

{ 

if((xarray[i] < x) ||(i=0)) 

{ 

y=-(pow( 1 0, 1 -((x-xarray [i+ 1 ])/(xarray [i]-xarray [i+ 1 ])))- 1 .0)/9.0*(yarTay[i]-yarray [i+ 1 D+yarray [i]; 
return y; 

}; 

}; 

}; 

return y; 

} 
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void funcinterpr(double yout[], double u[]) 
double xclose, exitpos; 


double y; 
double newx; 
double x; 

double* spiralexitx22,*spiralexitr22; 
int Ienspiralexit22; 

spiralexitx22=mxGetPr(mexGetGlobal("spiralexitx22")); 

spiralexitr22=mxGetPr(mexGetGlobal(”spiralexitr22")); 

lenspiralexit22=mxGetM(mexGetGlobal("spiralexitx22 M )); 

xclose=u[0]; 

exitpos=u[I]; 

x=u[2]; 

newx=x-exitpos; 
if (newx<0) 

{ 

y=1.0E+10; 

} 

else 

{ 

y-interpr(xclose-exitpos > spiralexitx 22 ,spiralexitr 22 ,lenspiralexit 22 )- 
if ((y>32766)&&(y<32768)) 

{ 

y=l .0E+10; 

}; 

}; 


yout[0]=y; 


} 

void mexFunction( 


int 

nlhs. 

Matrix *plhs[]. 


int 

nrhs. 

Matrix *prh$[] 



) 

double *yp; 
double *u; 
unsigned int m,n; 


/* Check for proper number of arguments */ 
if (nrhs !=1) { 

mexErrMsgTxt("Funcinterpr routine requires one vector input argument.")- 
} else if (nlhs > 1) { 

mexErrMsgTxt("Funcinterpr routine produces one vector output argument."); 
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/* Check the dimensions of u. */ 

m = mxGetM(UPIN); 
n = mxGetN(UP_IN); 

if (!mxIsNumeric(UP_IN) || mxIsComplex(UP_IN) || 
!mxIsFull(UP_IN) || !mx!sDouble(UP_IN) || 

(!(m == inputvec && n = 1) &&!(m = 1 && n == inputvec)) ) 

{ 

mexErrMsgTxtO'Funcinterpr routine requires 2 input vector."); 

} 


/* Create a Matrix for the return argument */ 

YP_OUT = mxCreateFull(l,l, REAL); 

/* Assign pointers to the various parameters */ 

yp = mxGetPr(YP_OUT); 
u = mxGetPr(UP_IN); 

/* Do the actual computations in a subroutine */ 

funcinterpr(yp,u); 

return; 
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FUNCNOSETIREFORCES.C is C function file which is converted to a CMEX file and 
called by ROTO SIMULINK picture FUNCNOSETIREFORCES in file 
MODELROTOBIG.M. It returns ground moment GM and nose gear forces FNOSE and 
MUNOSE. 

/* function [fhose,gm,MUNOSE] = funcnosetireforces(u)*/ 

/* function yout3 = funcnosetireforces(u)*/ 

#include <stdio.h> 

#include <math.h> 

#include "mex.h" 

/* Input Arguments */ 

#define UP_IN prhs[0] 

/* Output Arguments */ 

#define YP_OUT plhsfO] 

#define inputvec 13 


double fiincinterp(double x,double xarray[],double yaiTay[],int npts) 

double y; 
int i; 


y=-l-0; 

if (xarray[0] < xarray[npts-l]) 

{ 

for (i=l;i<=npts-l;i++) 

{ 

if «xarray[i] > x) ||(i==npts-l)) 

{ 

y = (x-xarray[i-l])/(xarray[i]-xarray[i-]])*(yarray[i]-yairay[i-l]) + yarray[i-l]; 
return y; 

}; 

}; 

} 

else 

{ 

for (i=npts-2;i>=0;i — ) 

{ 

if((xarray[i] > x)||(i=0)) 

{ 

y=(x-xarray[i+ 1 ])/(xarray [i]-xarray [i+ 1 ])*(yarray [i]-yarray[i+ 1 ])+yarray[i+ 1 ]; 
return y; 

>; 

}; 

} 

return y; 

} 

void funcnosetireforces(double yout[], double u[]) 
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int lenmuskidt; 

double U,sig,psin,pn,deldef,mumax,tvert,term2; 

double dell,delh,rl,rh,rl 1 l,r222,dell 1 l,de!222,hl 1 l,h222,terml,nl 1 l,n222; 
double phil 1 l,phi222,ql 1 l,q222,fhose,MUNOSE,deldel,stcs,tyaw,troll,ti!t,gm; 

double DT, SST ,K2 ,RB, OD ,SP ,HS, DELS, RP ,WS, S, MUROLL ,FLATMU ,D1. Kl. 
DELB; 

double tmpfy 111, tmpfy222, mulim 1 1st, mulim21st, fy 1 1 1 , fy222; 
double psisk2, vskid2 ,murat2 ,mulim2; 
double psiskl, vskidl, muratl ,muliml; 
double* MUSKIDT,*MUSKID; 

DT=*mxGetPr(mexGetGlobal("DT")); 

SST=*mxGetPr(TnexGetGlobal("SST M )); 

K2=*mxGetPi{mexGetGlobal("K2")); 

RB=*mxGetPr(mexGetGlobal( M RB n )); 

OD=*mxGetPr(mexGetGlobal("OD")); 

SP=*mxGetPr(mexGetGlobal("SP")); 

HS=*mxGetPr(mexGetGlobal( M HS n )); 

DELS=*mxGetPr(mexGetGlobal("DELS")); 

RP^mxGetP^mexGetGlobalC’RP")); 

WS=*mxGetPr(mexGetGlobal("WS")); 

S=*mxGetPr(mexGetGlobal("S")); 

MUROLL=*mxGetPr(mexGetGlobal("MUROLL")); 

FLATMU=*mxGetPr(mexGetGlobal("FLATMU M )); 

D 1 =*mxGetPr(mexGetGlobal( M D 1 ")); 

K 1 =*mxGetPr(mexGetGlobal("K 1 ")); 

DELB=*mxGetPr(rnexGetGlobal("DELB")); 

MUSKIDT=mxGetPr(mexGetGlobal( M MUSKIDT")); 

MUSKID^xGetPrOnexGetGIobalO'MUSKID")); 

lenmuskidt=mxGetM(mexGetGlobal( M MUSKIDT M )); 

/*not working, causes a memory problem 

tmpfy 1 1 l=*mxGetPr(mexGetGlobal("tmplyl 11”)); 
tmpfV222=*mxGetP^(mexGetGlobal( u tmp^V222 ,, )); 
mulim 1 lst=*mxGetPr(mexGetGlobal("mulirn 1 1st")); 
mulim21st=*mxGetPr(mexGetGlobal("mulim21st")); 
fy 1 1 1 =*mxGetPr(mexGetGlobal("fy 1 11")); 
fy222=*mxGetPr(mexGetGlobal("fy222")); 
mulim2=*mxGetPr(mexGetGlobal("mulim2")); 
mulim 1 =*mxGetPr(mexGetGlobal("mulim 1 "));*/ 

U=u[0]; 
sig=u[ 1 ]; 
mumax=u[2]; 
pn==u[3]; 
psin=u[4]; 

tmpfy 1 1 l=u[5]; 
tmpfy222=u[6]; 
mulim llst=u[7]; 
mulim21st=u[8]; 
fyl 1 l=u[9]; 
fy222*u[10]; 
mulim2=u[l 1]; 
mulim l=u[ 12]; 

deldef = fabs(SST*sin(sig/57.3)); 
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!*% initially set loads and deflections equal to zero */ 

/*% ('h' stands for high, T stands for low.)*/ 
dell = 0.0; 
delh = 0.0; 
rl = 0.0; 
rh = 0.0; 

/*% check whether or not the nose is off the ground. */ 
if(pn > 0.0) 

{ 

/*% assume all load is on one tire */ 
rh = pn; 
delh = rh/K2; 
if(rh >= RB) 

{ 

delh = rh/Kl+Dl; 

}; 

iffdelh > deldef) 

{ 

/*% assume both tires on high curve */ 
dell = D 1 -0.5 *(deldef-pn/K 1 ); 
rl = Kl*(dell-Dl); 
rh = pn-rl; 
delh = dell+deldef; 
iffdell < DELB) 

{ 

if(rh < RB) 

{ 

/*% both tires on low curve*/ 

dell = 0.5*(pn/K2-deldef); 
rl = K2*dell; 
delh = dell+deldef; 
rh - pn-rl; 

} 

else 

{ 

/*% one tire on low and one on high */ 

dell = (D 1 -deldef+pn/K 1 )/( 1 +K2/K 1 ); 
rl = dell*K2; 
delh = dell+deldef; 
rh = pn-rl; 

} 

} 

} 

}; 

/*% assign r and del according to steering direction*/ 
rl 1 1 = rl; 
r222 = rh; 
dell 11 = dell; 
del222 = delh; 
if(sig < 0.0) 

{ 

rl 1 1 = rh; 
r222 = rl; 
del 1 1 1 = delh; 
del222 = dell; 

}; 

h 1 1 1 = 0.85*OD*sqrt(del 1 1 1/OD - (del 1 1 l/OD)*(del 1 1 1/OD)); 

term 1 = SP+1 ,4/HS/OD*(dell 1 l*dell 1 1-DELS*DELS)*(SP+15.) + 0.44*RP; 
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if(dell 1 1/HS <= 0.338) 

{ 

nl 11 =terml*(WS*WS/3.32*(dell 1 1/HS- 1 .84*(del 1 1 l/HS)*(dell 1 1/HS))); 

> 

else 

{ 

n 1 1 1 = term 1 *(WS* WS/1 5.77*( 1 .-1 . 1 63*(del 1 1 1/HS))); 

}; 

h222 = 0.85*OD*sqrt(del222/OD - (del222/OD)*(del222/OD)); 

term2 = SP+1.4/HS/OD*(del222*del222-DELS*DELS)*(SP+15.) + 0.44*RP; 

if (de!222/HS <= 0.338) 

{ 

n222 =term2*(WS*WS/3.32*(del222/HS-1.84*(del222/HS)*(del222/HS))); 

} 

else 

{ 

n222 = term2*(WS*WS/15.77*(l.-1.163*(del222/HS))); 

}; 

if(n 1 1 1 > 0.0) 

{ 

psiskl = psin - fyl 1 1/nl 11; 

vskidl = fabs(U/1.689*(sin(psiskl/57.3))); 

muratl = funcinterp(vskidl,MUSKIDT,MUSKID,lenmuskidt); 

if (FLATMU == 1) 

{ 

muratl = 1.0; 

}; 

mulim 1 = mumax*muratl; 

} 

if(n222 > 0.0) 

{ 

psisk2 = psin - fy222/n222; 

vskid2 = fabs(U/1.689*(sin(psisk2/57.3))); 

murat2 = funcinterp(vskid2,MUSKlDT,MUSKID,lenmuskidt); 

if (FLATMU == 1) 

{ 

murat2 = 1.0; 

}; 

mulim2 = mumax*murat2; 

}; 

if (rill >0.0) 

{ 

phil 11 = nl 1 l*fabs(psin)/muliml/rl 11; 
if(phi 1 1 1 <= 0.1) 

{ 

ql 1 1 = 0.8*hl 1 1/(1 .0-. 1482*phi 1 1 l*phil 1 1); 

} 

else 

{ 

if (phil 11 <= 0.55) 

{ 

ql 1 1 = hi 1 l*(phil 1 1-phil 1 l*phil 1 1 -0.0 1 )/(phi 1 1 l-0.1482*phil 1 l*phil 1 l*phil 11); 

} 

else 

{ 

if (phi 111 <= 1.5) 

{ 

ql 1 1 = hi 1 1 *(0.2925-0. l’phil 1 l)/(phil 1 l-0.1482*phil 1 l*phil 1 l*phil 1 1); 
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} 

else 

{ 

qlll = hi 1 1 *(0.2925-0. l*phil 1 1); 

}; 

} 

} 

} 

else 

{ 

phi 111 =0.0; 
qlll = 0.0; 

}; 

fyl 1 1 =muliml*rl 1 l*(phil 1 l-0.1482*phil 1 l*phil 1 l*phil 11); 
if(phil 1 1 > 1.5) 

{ 

fyl 1 1 = muliml*rl 1 1; 

}; 

if(psin < 0.0) 

{ 

fyl 1 1 = -fyl 11; 

}; 

if(r222 > 0.0) 

phi222 = n222*fabs(psin)/mulim2/r222; 
if(phi222 <= 0.1) 

{ 

q222 = 0.8*h222/(1.0-.1482*phi222*phi222); 

} 

else 

{ 

if (phi222 <= 0.55) 

{ 

q222 = h222*(phi222-phi222*phi222-0.01)/(phi222-0.1482*phi222*phi222*phi222); 

else 

{ 

if(phi222 <= 1.5) 

{ 

^q222 = h222 *(0.2925-0. 1 *phi222)/(phi222-0. 1482*phi222*phi222*phi222); 

else 

{ 

q222 = h222*(0.2925-0. 1 *phi222); 

}; 

} 

} 

} 

else 

{ 

phi222 = 0.0; 
q222 = 0.0; 

}; 

fy222 =mulim2*r222*(phi222-0. 1 482*phi222*phi222*phi222); 
if (phi222 > 1 .5) 

{ 

fy222 =mulim2*r222; 

}; 
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if(psin < 0.0) 

{ 

fy222 = -fy222; 

}; 

if(U > 0) 

{ 

/*% relaxation lengtHS if mumax increasing*/ 
if(mulim 1 * 1 .05>=mulim 1 1st) 

{ 

tmpfy 1 1 1 =fy 1 1 1 +( tmpfy 1 1 1-fyl 1 l)*exp(-DT/(3/U)); 
fyl 1 l=tmpfyl 1 1; 

} 

else 

{ 

tmpfy 1 1 1 =fy 111; 

}; 

mulimllst=muliml; 
if(mulim2* 1 ,05>=mulim21st) 

{ 

tmpfy222=fy222+(tmpfy222-fy222)*exp(-DT/(3/U)); 

fy222=tmpfy222; 

} 

else 

{ 

tmpfy 222=fy222; 

}; 

mulim21st=rnulim2; 

}; 

fhose=fy 1 1 1 +fy222; 

MUNOSE=fnose/pn; 

/*% GROUND MOMENT ON NOSEWHEEL STRUT */ 
deldel = SST*sin(sig/57.3); 
stcs = SST/S*cos(sig/57.3); 

tyaw = fyl 1 l*((OD/2.-dell 1 l)*stcs + ql 1 1) + fy222*((OD/2.-del222)*stcs + q222); 
troll = MUROLL*(rl 1 l-r222)*S; 
tilt = asin(deldel/S); 

tvert = r222 * stcs* (0 . 5 * S-(0.5* OD-del222)*tan(ti lt))-r 1 1 l*stcs*(0.5*S+(0.5*OD-dell 1 l)*tan(tilt)); 
gm = tyaw + troll + tvert; 

/*for some reason this is not working, causing memory problems 

mxSetPr(mexGetGlobal("tmpfyl Il"),&tmpfyl 1 1); 
mxSetPrfm exGetG lobal("tmpfy222 ”),&tmpfy 222); 
mxSetPr(mexGetGlobal("mulimllst"),&mulimllst); 
mxSetPr(mexGetGlobal("mulim21st"),&mulim21st); 
mxSetPr(mexGetGlobal("fyl 1 l"),&fyl 1 1); 
mxSetPr(mexGetGlobal(”fy222"),&f\222); 
mxSetPr(mexGetGlobal("mulim2' , ),&mulim2); 
mxSetPr(mexGetGlobal("muliml"),&muliml);*/ 


yout[0]=fnose; 
yout[l]=gm; 
yout[2]=MUNOSE; 
yout[3]=tmpfyl 11; 
yout[4]=tmpfy222; 
yout[5]=mulimllst; 
yout[6]=mulim21st; 
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yout[7]=fyl 11; 
yout[8]=fy222; 
yout[9]=mulim2; 
yout[10]=muliml; 


} 


void mexFunction( 
int 

Matrix *plhs[], 
int 

Matrix *prhs[] 

) 

{ 

double *yp; 
double *u; 
unsigned int 


nlhs, 

nrhs, 


m,n; 


/* Check for proper number of arguments */ 
if (nrhs !=]) { 

mexErrMsgTxt("Funcnosetireforces routine requires one vector input argument "V 
} else if (nlhs > 1) { 

mexErrMsgTxt( 'Funcnosetireforces routine produces one vector output argument."); 


/* Check the dimensions of u. */ 

m = mxGetM(UPIN); 
n = mxGetN(UPIN); 

if ( ImxIsN umeric(U P_fN) || mxIsComplex(UP_IN) || 

!mxIsFull(UP_IN) || !mxlsDouble(UP_IN) || 

(!(m = inputvec && n == 1) &&!(m = 1 && n = inputvec))) 

mexErrMsgT xt("F uncnosetireforces routine requires a 13 input vector."); 


/* Create a Matrix for the return argument */ 

YP OUT = mxCreateFull(l,l 1, REAL); 

/* Assign pointers to the various parameters */ 

yp = mxGetPr(YP_OUT); 
u = mxGetPr(UP_IN); 

/* Do the actual computations in a subroutine */ 

funcnosetireforces(yp,u); 

return; 


A- 62 



MODELROTO 



A-63 


NEWNAVDATfr / EXITPOS 








AERO DPS 



A-65 

























A-66 


G & Jerk 





















A— 68 




A— 69 






A-70 


















A-71 


►M DXN [oJJ Relational I — nFk - 

Muxl Constant Operator algebraic loop 























YAW ANGLE OF L&R MA IN GEAR i i r»|u[3]-u[1]/57,3*sqrt(pow((TREAD/2),2)+pow(B,2))*sin(atan(TREAD/2/B)-uf2V57.3) 
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To Workspace 1 
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CTRL LAWS 


see RUNROTO for inputs to predict exit selection logic 


lateral displacement, yioc 

~5~] laL dis P ratp V* nr 
yaw rate HT 

centerline r adius 

Iruddef) 

Constant, rudder at touch down 

/newnavdata, 
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/ udropout 




navigation al X p os., xnav 

S— — 

navigational V pos., ynav 


longitudinal gmd speed, UNAV 
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eg distance require d to dear runway, clrdist 
ASYMBRAKEl 


Constant 

n_ 

elapsed runway occupancy time, runtime 


main gear brake press., 
~ i, PBC, PfeL 
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thrust Ctrl law 


autoreverse 
thrust command 
multiplier 
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APPENDIX A: EVENT TIMELINE FOR A REPRESENTATIVE ROTO IMPLEMENTATION - 
MANUAL LANDING / AUTOMATIC ROTO DECELERATION AND RUNWAY EXIT 
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EVE ^!n E V , . N ^™“ * “ERfS^TATIVE ROTO IMPLEMENTATION (CONTINUED): 
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EVENT TIMELINE FOR A REPRESENTATIVE ROTO IMPLEMENTATION (CONTINUED): 
MANUAL LANDING / AUTOMATIC ROTO DECELERATION AND RUNWAY EXIT 
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EVENT TIMELINE FOR A REPRESENTATIVE ROTO IMPLEMENTATION (CONTINUED): 
MANUAL LANDING / AUTOMATIC ROTO DECELERATION AND RUNWAY EXIT 
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EVENT TIMELINE FOR A REPRESENTATIVE ROTO IMPLEMENTATION (CONTINUED): 
MANUAL LANDING / AUTOMATIC ROTO DECELERATION AND RUNWAY EXIT 
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